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1. Introduction. The purpose of our lecture is, briefly stated, the discus- 
sion of periodic solutions of the differential equation 


# + f(x) = F cos ot, (1) 


in which f(x) is in general a nonlinear function of x. For f(x) =g/l sin x we 
have, for example, the case of the simple pendulum with a periodic external 
force. Such a differential equation may occur in electrical problems when iron 
is contained in the magnetic circuit.! Any system with one degree of freedom? 
and consisting of a mass and an elastic restoring force will lead to a differ- 
ential equation of the type (1) if the displacement x is not kept very small. 

The problems to be discussed here differ essentially from those of self- 
excited oscillations where negative damping is involved, as is the case for the 
group of problems first treated by van der Pol. Neither the nonlinear char- 
acter of the restoring force as such nor linearity or nonlinearity in possible 
friction forces is responsible for the differences. In the case of self-excited 
oscillations results are altered but little if a nonlinear restoring force is as- 
sumed, while the results to be discussed here would not be altered in general 
character if nonlinear positive damping were assumed. In other words, it ap- 
pears that even a slight negative damping will dominate the entire phe- 
nomena. 

Our discussion will yield, for the most part, results which are already well 
known. Much of the material can be found in the standard textbooks, e.g., 
in the books of Timoshenko [28] and den Hartog [5] cf. also [11, 20]. How- 


* Lecture given at the Conference on Non-Linear Mechanics, Brown University, August 
1-3, 1942. Manuscript received Oct. 15, 1942. 

1 For this subject see [32] and the list of references given there. 

2 Some cases of free oscillations with two degrees of freedom have been treated [4] and 
[12] (3. Mitt.). For experimental results in a case with infinitely many degrees of freedom, 
see [22]. 
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ever, the problems are of sufficient intrinsic importance that they merit 
thorough consideration and discussion. 

We are interested, to begin with, in obtaining insight into the physical 
phenomena connected with-the case under discussion. It turns out that it is 
highly useful to begin by obtaining information of a qualitative character. 
Once this has been done, it is not very difficult to see how one should proceed 
to obtain accurate quantitative information. The validity of the methods to 
be discussed here could be justified rigorously in most cases, but we donot 
wish to take up these matters here. We prefer to stress methods which lead 
to the essential qualitative information in as simple a manner as possible. 
What might be found original in this lecture consists largely in the exploita- 
tion of this point of view. 

2. Nonlinear restoring forces. The differential equation (1) has been 
treated for a considerable number of different restoring forces f(x), cf. [5, 6. 7. 
16, 24, 25, 28]. Some of these cases are indicated in fig. 1: 
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The forces indicated in fig. 1 are symmetrical, i.e. f(—x) = —f(x). This need 


not always be the case, of course; but we shall assume it in what follows. The 
last two examples, for the hard and soft springs, indicate the essential distinc- 
tion between different types of springs, i.e., those for which the stiffness f'(x) 
is an increasing or a decreasing function of x. The first case is typified by a 
mass attached to a stretched string, the second by the pendulum. 
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Since methods and also qualitative results appear not to depend greatly 
upon the special form of f(x), we shall choose for this function always the 
expression 


f(x) = ax + Bz’, a> 0. 


In any case there is no great loss in generality involved for moderate values 
of x, since our expression might be regarded as the linear and cubic terms in 
the power series for f(x). Note also that 8>0 characterizes a hard spring 
(stretched string), while 8<0 characterizes a soft spring (pendulum). 

The differential equation (first treated by Duffing [8]) which we want to 
consider in the following is, thus, 


% + (ax + Bx*) = F cos at. (2) 


3. Periodic solution. We seek periodic solutions of (2) with frequency w. 
(The term frequency will be used here in place of the more correct term circular 
frequency, since no confusion is likely to result). That solutions of (2) other 
than periodic ones exist is certainly true (even unbounded ones exist for 8 <0). 
However, the literature of the subject is almost entirely devoted to the peri- 
odic solutions, which are the interesting ones in the engineering applications. 
Apparently the experimenters always find periodic motions, at least after 
some transient motions have died out. Viscous damping, which is always 
present in any actual case, seems to act in such a way that the motions in a 
wide variety of cases tend, as > ©, to periodic ones. It would be of consider- 
able interest to prove that the solutions of (2) with a damping term added 
(and under appropriate conditions) are of this character. 

We turn, then, to the problem of finding periodic solutions of (2). No ex- 
plicit solutions of this differential equation are known and we are forced to 
turn to approximate methods. Perhaps the simplest of these is the iteration 
method. Let us write (2) in the form 


= — (ax + Bx*) + F cos wl (3) 
and insert 
xo = A cos wt (4) 


as a first approximation in the right-hand side. This means, in effect, that we 
assume 8 to be small so that a motion not greatly different from a simple 
harmonic motion can be expected. In fact, we assume in all of our methods 
but one that @ is small. Upon using the identity 


cos* wt = 3 cos wt + } cos 3wt (5) 
we find 
#, = — (aA + 38A* — F) cos wt — 1 8A* cos 3wt (6) 


as equation for the next approximation x. Integrating this one obtains 
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BA 
= cos 3wt, (7) 


w" 





1 1 
x = — (ad + 364? — F) cos wt + — 
w? 36 


the integration constants being taken as zero to insure that x; and the next 
following approximation x2 be periodic. 

So far this is quite straightforward. What to do from now on is not so 
clear. One might seek to continue the process by inserting x; in the right side 
of (3) in order to find an x2, etc., which would be a rather natural procedure. 
The significance of this procedure can be illustrated by the following figure 
which indicates the well known response curve for the linear forced vibration, 
i.e., a curve showing the amplitude | A | of the steady forced vibration asa func- 
tion of the frequency w, the force amplitude F being a parameter (fig. 2): 
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The procedure just outlined would mean that w is held fixed while | A| 
is left open, for 80. This would yield certain curve (shown dotted) in the 
neighborhood of the linear response curve. It is clear that there would be 
difficulties near w=+/a. But worst of all, the really essential features of the 
response curves for 8#0 would not be obtained at all. 

The method used by Duffing [8], who seems to have been the first to 
obtain the significant results for the differential equation (2) was as follows: 
the coefficient A, of cos wt in (7) is taken equal to A in (4), on the ground that 
A; should differ but little from A if (4) is truly a reasonable first approxima- 
tion. Also, Duffing argues, such a procedure would furnish the exact result 
in the linear case (@=0) and might hence be expected to yield good results 
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for 8 small, which we assume here anyway so that convergence can be ex- 
pected. This reasoning of Duffing leads thus to what proves to be the decisive 


relation (for A #0) 
F 
wt = a + BAP — (8) 


between the amplitude A and frequency w of the periodic solution. The rela- 
tion (8) has been written purposely so that w* is given as a function of A: as 
we shall see it is decisive to consider A as a prescribed quantity in terms of 
which w is to be determined. 

Before discussing the significance and interpretation of (8), let us first 
obtain the same relation by a procedure which is perhaps a more systematic 
one. Equation (3) is modified by adding a term w*x to both sides: 


t+ wx = — [(a — w*)x + Bx*] + F cos ot. (9) 

As a first approximation x9 to a periodic solution we begin with the solution 

of (9) for a=w?, 8=0, F=0, i.e., with a free undamped linear oscillation. 
This leads at once to 

xo = A cos wt, (10) 


with A arbitrary. Upon insertion in the right-hand side of (9) and use of (5) 


once more we find 
1 + wx, = {(w? — a)A — 38A*+F} cos wt — 48A* cos 3wt. = (11) 


We require always that the solution be periodic; hence it is necessary that the 
coefficient P; of cos wf in the right hand side of (11) should vanish in order to 
avoid the resonance case and hence the occurrence of nonperiodic terms in 
the solution of (11). The vanishing of this coefficient yields, obviously, the 
same relation as (8), which is to be regarded as an equation to determine w 
after A has been prescribed. Once this relation has been satisfied, the solution 
of (11) will be 


cos 3wi, (12) 





pA 
x; = A, cos wt + B, sin wi + 


32w? 
in which A; and B, are arbitrary. The method of fixing A; and B, is, however, 
now clear. We simply set A1=A and B,=0, the value of w being left open. 
If we stop with the degree of approximation implied in x; we should have, 
then, as approximate solution 


1 
x = A cos wt + te cos 3wt. (13) 
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It is perhaps worth while to consider how the iterations should proceed: 
x1 from (13) is inserted again in the right-hand side of (9). A term P2 cos wt will 
again occur and we must require, as before, that P2=0. The equation P;=0 
yields an improved relation between w and A that takes the place of (8). 
After integration the coefficient of cos wt is again prescribed to have the value 
A while the coefficient of sin wt is taken as zero. 

One could interpret this method as meaning that the first coefficients in 
the Fourier series for x(t) are prescribed once and for all to be A and zero, 
while the frequency is left open. In what follows it is convenient to speak of 
the quantity A as the amplitude of the oscillation, though it should be referred 
to rather as the first Fourier coefficient. 

Our procedures, in any case, require that A be fixed while w, the frequency, 
is to be determined by relation (8). Why this at first sight seemingly unnatu- 
ral procedure leads to the desired results can be best seen through a discussion 
of relation (8), which yields a set of curves, the response curves (in first ap- 
proximation), with the force amplitude F as parameter. Usually the sign of 
A is not considered essential (it means a phase shift) so that only |A| is 
plotted. The curves for | A| against w are then readily seen to be of the form 
shown in fig. 3, where curves for 8>0, 8=0 and 6 <0 are given: 
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The interpretation of the effect of the nonlinearity on the response curves is 
quite clear: The entire family of response curves (with F as parameter) for the 
linear forced oscillations is bent left or right depending upon whether the 
spring is a soft or a hard one. That these curves mirror the essential facts, 
at least for 8 and A not too large, cannot be doubted. 

We can now see the reasons why A, instead of w, should be prescribed. To 
begin with, we started our approximation with w=/a, B=0, F=0, that is, 
at point 1 in fig. 4, which superimposes curves for 8=0 upon curves for 
8>0. One sees that, for 8B=0, A is completely arbitrary while w is fixed 
and it is therefore necessary to prescribe A but leave w open for B¥0. 
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The further steps in the iterations correspond to the passage from this 
state to one indicated by point 2 in fig. 4. One sees also that to hold w fixed 
while varying A could yield at most one of the branches of the response 
curves. 

This discussion indicates that the guiding principle to be followed inde- 
pendent of the particular scheme of approximation used is to observe what is 
arbitrary and what is fixed in the linear case, i.e., for 8B=0. For 80, any 
quantity that is arbitrary for 8=0 should be prescribed, and vice versa. 
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4. Other methods. Before proceeding to draw the interesting physical con- 
clusions from the response curves, we prefer to discuss some other methods 
of approximating the periodic solutions of (2). 

One of the best known and most often used methods for nonlinear prob- 
lems is the perturbation method, cf. [28], which consists for our problems 
in developing the solution x(t) in powers of a parameter, say €, with coeffi- 
cients prescribed to be periodic functions of ¢: 


x= xX texx t+ Cext+-:-- P (14) 


However, as we have seen, the frequency should not be prescribed but rather 
left open, and this requires that w should also be developed with respect to e: 


w= wo + ew + wet ---, (15) 


where the w; are constants. It is awkward to work with functions which have 
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a variable period, but one can avoid this by introducing a new variable 0 =wt, 


so that (2) becomes 


wx” + (ax + Bx*) = F cos 8, (16) 


the prime meaning differentiation with respect to 6. The functions x; in (14) 
can now be prescribed to have the fixed period 27. The perturbation series 
(14) and (15) are now to be inserted in (16), and the coefficients of the powers 
of € in the resulting series are equated to zero. The result is a set of linear 
differential equations for the x;. The parameter € is to a certain degree arbi- 
trary. A very reasonable choice for it would be to take e=8, so that (14) could 
be considered a development in the neighborhood of the solutions of the 
linearized vibration problem. The first few steps lead to essentially the same 
result as those obtained above by iteration. One finds, for example, that rela- 
tion (8) results from the periodicity requirement. It should be said that the 
convergence of the perturbation series could be established with no great 
difficulty. 

Since the solutions which interest us are periodic ones, they will possess 
Fourier series developments. This suggests writing the solutions x(t) at the 
outset as such a series with undetermined coefficients: 


x = Ap +A, coswt + Azcos2wi+-:-:- 
+ B, sin wt + Basin 2wi+--:. 


(17) 


Upon insertion into (2) one obtains an infinite set of nonlinear algebraic equa- 
tions for the A; and B; which in their turn are to be solved by iterations or 
perturbations. Our previous experience indicates, however, that we should 
not prescribe w but rather the first coefficients in (17) and this is decisive in 
fixing the correct set-up for iterations, say. Probably the Fourier series method 
is the simplest for actual computation once one has the clue to the correct 
procedure. 

The approximation methods considered so far have all taken the solution 
for the linear free oscillation as a starting point. Although it is true that the 
significant qualitative features of the response curves can thus be obtained, 
it is clear that other procedures could be expected to yield more accurate 
results. Instead of beginning with the free linear oscillation, for example, one 
might begin with the free nonlinear oscillation, that is, at point 1’ rather 
than point 1 in fig. 4. For a given value of the force amplitude F it is 
quite obvious that the solution corresponding to point 1’ is much nearer to 
the desired solution (point 2) than is point 1. It is, of course, quite feasible 
to begin with a solution corresponding to point 1’, i.e., with a solution of 
#+ f(x) =0, since such a differential equation can be solved by explicit in- 
tegration. With this solution as a basis one might then proceed by the per- 
turbation method using F as small parameter, as for example has been done 
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by G. Hamel [9] and R. Iglisch [12, 13]. The method of variation of parame- 
ters has also been used, for example, by E. W. Brown [4], and A. N, Krylov 
and N. Bogolitiboff [18, 19] have used such methods in a considerable number 
of papers. 

The method of M. Rauscher [24, 25] is an iteration procedure which also 
takes the solution corresponding to point 1’ as its starting point. It is con- 
venient here also to introduce the variable 6 =w# in order to work with func- 
tions which have a fixed period (as in the perturbation method). With this 
variable our differential equation is 


wx” + f(x) = F cos 0. (18) 


The idea of Rauscher is to take as first approximation the periodic solution 
xo(0) of 
wx” + f(x) = 0 (19) 


which has a prescribed amplitude A (i.e., (Xo)max =A). The quantity wo is to 
be determined in such a way that xp will have the period 27. This is a problem 
that can be solved by explicit integration, which may be carried out graphi- 
cally if necessary. The quantity wo is, of course, taken as a first approximation 
to the unknown frequency w. The essentially new idea of Rauscher consists 
then in replacing @ in the right hand side of (18) by the function 69(x) ob- 
tained by inverting the solution x9(@) of (19), over a half period in which 
xo is a monotone function of @. This yields as equation for the next approxima- 


tion: 


wix"” + f(x) — F cos [6(x)] = 0, (20) 


which is an equation of the same type as (19): it differs only in that the re- 
storing force is different. Its solution, for the same fixed A and same period 27, 
yields a second approximation x;(@) and a new value @ for w, which will in 
general be different from wo. The general scheme of this iteration method thus 
consists in finding the solution x,(0@) of period 27 of 


wrx" + f(x) — F cos [6,-:(x)] = 0 (21) 


for (%n)max =A, where 6,_:1(x) has been obtained by inverting the solution of 
the equation for (n—1). At the same time the mth approximation a, for is 
found. Actually, there is no need to invert the solution each time: The 
Rauscher method makes direct use of the fact that the solution of an equation 
of the type (19) is obtained most conveniently at the outset in the form 
6 =0(x). In the case of the pendulum, for example, the time is obtained as an 
elliptic integral which must be inverted to obtain the displacement as an 


elliptic function of the time. 
The amount of numerical computation is not small in this method, but 
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the convergence is so rapid that it may well be the best method of obtaining 
accurate quantitative results. Rauscher finds that the second step usually 
furnishes a very close approximation. It might be objected that the method 
of Rauscher could not be used if viscous damping were present since the differ- 


ential equation 

x” + cx’ + f(x) =0 
cannot be integrated explicitly. However, the term cx’ can be treated in 
much the same manner as the term F cos @ once a phase shift has been in- 


troduced. 
Finally, brief mention should be made of two other methods. The Ritz 


method has been used [8, 9, 27, 31 ]. It consists in minimizing the expression 


r B 
f 1** — = x* + 2xF cos ot a 
P 2 


with respect to periodic functions x(/) of period T =27/w. The expression at- 
tains a minimum, however, only for soft springs (8 <0). 
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In case one is interested in obtaining a solution satisfying given initial 
conditions (such a solution need not be periodic), the graphical method of 
E. Meissner [23] is to be recommended. The method is particularly well de- 
signed for just the type of differential equation considered here. 

5. Jump phenomena. We turn now to the physical conclusions which can 
be drawn from the form of the response curves. Instead of discussing these 
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phenomena on the basis of the curves of fig. 3, we prefer to consider the 
analogous curves which result under the assumption of slight viscous damp- 
ing in the system. By analogy with the response curves in the linear case with 
damping, one would expect the curves to be rounded off as indicated in fig. 5: 
This is a point to which we shall return a little later. 

In the following fig. 6, we show a typical response curve (i.e., for a given 
value of the force amplitude F) for the case of a hard and a soft spring: 
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Let us imagine an experiment performed in which the amplitude F of 
the external force is held constant, but its frequency w is slowly varied. (If 
w were held fixed and F slowly varied, the results would be of the same gen- 
eral character.) Consider first the case of a hard spring (8 >0) and suppose 
that w is rather large at the beginning of our experiment, i.e., we start at 
point 1 on the curve. As w is decreased A slowly increases through point 2 
until point 3 is reached. Since F is held constant, a further decrease in w 
would require a jump from point 3 to point 5 with an accompanying increase 
in the amplitude A, after which A decreases with w. Upon performing the 
experiment in the other direction, i.e., starting at point 6 and increasing w, 
the amplitude follows the 6-5-4 portion of the curve, then jumps to point 2 
and afterwards slowly decreases. The circumstances are quite similar with a 
soft spring, but the jumps in amplitude take place in the reverse direction. 
It would not have been necessary, we see now, to consider the influence of 
damping in order to conclude that a jump from point 3 to point 5 should 
take place (for 8>0) on decreasing w, but the jump from point 4 to point 2 
on increasing w would be inexplicable on the basis of the curves of fig. 3. 

This curious behavior of nonlinear systems has often been observed by 
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experiment. Probably the first to discuss it experimentally was O. Martiens- 
sen [21]. Martienssen carried out his experiments with an electrical apparatus 
involving a condenser and an inductance with an iron core; he ascribed the 
results to the nonlinearity. 

6. Effect of viscous damping. We wish to discuss briefly a method of 
treating the effect of viscous damping, i.e., the modifications which result upon 
adding a term cx, c>0, to, the left-hand side of (3): 


+ cé + (ax + Bx’) = F cos wt —G sin at. (22) 


The amplitude H=1/ F?+G of the external force is prescribed, but the ratio 
F/G is not prescribed in order to permit a phase difference between the ap- 
plied force and the resulting periodic displacement. 

The iteration method can be applied in much the same manner as above 
by beginning with x»=A cos wf as a first approximation. In place of the 
relation (8), we now have the following two relations: 


(a — w*)A + 38A* = F, (23) 
cwA =G. (24) 


Once can draw an interesting, though not unexpected, conclusion immedi- 
ately. It is that no periodic motion except the state of rest can result if there 
is damping but no external force, i.e., if F=G=0 with c#0. In fact, (24) is 
satisfied in this case only for A =0 which implies x»=0. 

Equations (23) and (24) are squared and added to obtain 


[(a — w*)A + 38A3]? + c2w?A? = F2 4G? = H?, (25) 


which gives the relation between the amplitude A and frequency w of the 
vibration which results from the force amplitude 7. Equation (25) replaces 
(8) for this case; it of course reduces to (8) when c=0. 

The discussion of the response curves furnished by (25) turns out to be 
quite easy and to yield the expected results. One finds that the curves have 
only one branch instead of two as in the undamped case, and that each curve 
shows a single maximum for A. The locus of the maxima (for different H) is a 
curve which runs close to the response curve for the free undamped oscillation 
(H=c=0), always remaining to the left of it. With no damping there is only 
one vertical tangent on each response curve. With damping a second vertical 
tangent appears; the new point of tangency lies near to the curve for free 
undamped oscillations. These results have already been indicated in figs. 5 
and 6. 

It is of some interest to note that the curves for values of H less than a 
certain limit value will have no vertical tangents, which means that for such 
amplitudes no jump phenomena would occur. 

The effect of viscous damping on the response curves appears to have been 
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discussed first by E. V. Appleton [1]. Appleton uses the method of variation 
of parameters. He investigates also in detail the stability of the solutions. It is 
perhaps of interest to note that Appleton was led to the problem through ob- 
serving the peculiar behavior of a certain galvanometer at the Cavendish 
laboratory at Cambridge. This behavior corresponds to what was called jump 
phenomena above. 

7. Stability questions. We shall discuss briefly the stability question for 
the cases without damping. Probably the most difficult question in any con- 
sideration of dynamic stability is that of the definition of stability itself in a 
reasonable way. Into this question we do not wish to enter here. We take, 
simply, the following often used definition of stability: let x(¢) and x(t) + 6x(é) 
be two solutions of our differential equation. Consider the variational equation 
which results when x+ 6x is inserted in the original differential equation and 
powers of 46x above the first are neglected. If all solutions 6x of this equation 
are bounded, then x(¢) is said to be stable, otherwise unstable. 

In particular, the variational equation for Duffing’s equation (2) is 


5% + (a + 38x*)ix = 0. (26) 


In the cases under discussion here x(¢) is a periodic function and thus (26) 
is a Hill’s equation. The question is then whether for a given x(t) all solutions 
of (26) are bounded or not. Or we might put the problem a little differently 
by referring to the A-w-plane, the points of which characterize the periodic 
solutions x(#). The problem thus is to divide this plane into regions which 
correspond to stable solutions and others which correspond to unstable ones. 
Questions of this kind have been discussed, cf. for example [1, 18, 29, 31], 
but the problem has apparently not been solved completely. 

Some conclusions can, however, be reached without much difficulty. Con- 
sider, for example, fig. 7, which indicates the response curve for F=0 (i.e., 
for the free oscillation) together with the locus of vertical tangents on the 
response curves for F>0. (In connection with the latter curve see [12, 14, 
15]). Our discussion of the jump phenomena would lead us to expect that the 
shaded area between these two curves corresponds to solutions which are 
unstable in some sense—quite possibly in the sense of the above definition. 
It can at least be shown that the two curves in question really do represent 
boundary curves between regions of stable and unstable solutions in the sense 





of the definition. 
If one were to be contented with taking A cos wt as a sufficiently close 


approximation to x(t), then equation (26) would be a Mathieu equation. 
For the special case of the Mathieu equation the problem of separating stable 
and unstable solutions has been completely solved. If one makes use of this 
theory it turns out that our two curves have contact of the second order at 
w=/a, |A| =0 with curves delimiting a region of unstable solutions of the 
Mathieu equation, the unstable region being that within the cusp. 
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8. Subharmonic response. Up to now we have considered periodic solu- 
tions of (3) for which the frequency is always exactly the same as that 
of the external force F cos wt. Permanent oscillations with a frequency 
3,4,---1/n, +--+ of that of the applied force can, however, occur in non- 
linear systems, in particular in our case of the Duffing equation. To this phe- 
nomenon the term subharmonic resonance is usually applied.* 
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The fact that subharmonic oscillations occur can hardly be denied since 
they have been so often observed. But it is not an entirely simple matter to 
give a plausible physical explanation for their occurrence. Let us recall the 
behavior of linear systems. If the frequency of the free oscillation of a linear 
system is w/n (n an integer, say) then a periodic external force of frequency w 
can excite the free oscillation in addition to the forced oscillation of fre- 
quency w. But since some damping is always present the free oscillation is 
damped out so that the eventual steady state consists solely of the oscillation 
of frequency w. Why should the situation be different in a nonlinear system? 
The explanation usually offered is as follows: Any oscillation of a nonlinear 
system contains the higher harmonics in profusion. Hence it is possible that 
an external force with a frequency the same as one of these might be able to 
excite and sustain the harmonic of lowest frequency. Of course this requires 
that the damping (more precisely the ratio c/H) be not too great and that 
proper precautions of various kinds be taken. 

We shall not attempt to present a solution of the problem of subharmonic 


3 For literature, see [2, 3, 17, 18]. 
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response for the Duffing equation in all generality. Rather, we shall treat only 
one special case. Also, damping will be neglected even though it should be 
shown that subharmonic response occurs in spite of damping. This could be 
done, but it would complicate our calculations unnecessarily without causing 
any significant qualitative changes in the results, particularly if the damping 
is small. In any case we are interested here, as in our previous discussion, in 
qualitative rather than accurate quantitative results. 

There is some advantage in introducing once more the variable 6 =wt as 
the new independent variable. The differential equation is 


wx” + (ax + Bx*) = F cos 8, (27) 


in which it is to be remembered that w represents the frequency of the applied 
force F cos wt. Our object is to find a periodic solution with frequency w/n. 
We restrict ourselves to the case »=3. We apply the Fourier series method 
and set 


6 50 
x = Ays —s + Aol + Aon + tee 


by analogy with the usual or harmonic case. (Terms in even multiples of 6/3 
and sine terms drop out). Upon substitution of the series in (2) one finds the 
relations 


w? 3 2 aS 
(« Pe =)A v3 + 38[A 3 + AiyjsAi + 2Arj3Ai + -- J at (28) 


(a — w )Ay + 38lAvs + GAyeAr + 34:+---) =F. (29) 


Equations (28) and (29) take the place of equation (8) which was fundamental 
for the harmonic case. 

Again we are faced with a problem of interpretation. Our previous experi- 
ence offers a guide: we set B = 0 (i.e., we consider the linear case) and observe 
that A; must be taken zero unless w=3+/a. If, however, w=3+/a then Ays 
can be taken arbitrarily, while A; = —(1/8a&) F is determined for a given value 
of F. The term Aj;3 cos 6/3 is evidently the free oscillation of arbitrary am- 
plitude which may be superimposed on the forced vibration —(F/8w) cos @ 
in the linear case. Hence we should prescribe A1/3 arbitrarily for 840 but 
then hold it fixed. The quantity F should also be prescribed, but the quanti- 
ties w and A, which were fixed for 8 =0 should now be considered as functions 
of Ay3 and F. 

With this in mind we rewrite equations (28) and (29) in slightly different 
form in order to perform iterations conveniently. This yields, after a division 
of (28) by A130: 
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2 . 2 a . 

w = 9a + 228[Ays + Arjsdi + 241+ --> |, (30) 
, 2 3 2 3 

— 8aA,; =F + (w — 9a)A; — 38[A1s + 641341 + 34, +--- |, 


or 


; 3 2 Ps oo 78 ‘ 
— 8aA; = F — 38[Ays — 21A1241 — 27AjzA1 — 51A, +--+: |], © (31) 
the last equation resulting through elimination of w?. 
We begin the iterations with the values for 8 =0, i.e., with A1,3 prescribed, 
w =3/a, and A;= — F/8a=A. The next step yields 


2 2 
w = 9a + 228[A13 + AyysA + 2A |, (32) 
A+ Se Miho ~ HbA ~ 541, (33) 


A, 


Equation (32) is readily discussed: it represents an ellipse or a hyperbola 
in an w-A}4/3-plane, depending on the sign of 8. Also w has an extremum for 
Aij3= —34A and w? has as value there 


w* = 9(a + #8A?). (34) 


Thus the subharmonic vibration can exist only if 





w $ 3Va+FhBA? for (35) 


If 80, we may conclude that no vibration with w=3+/a can exist; i.e., no 
subharmonic response with exactly the frequency of the linearized problem 
can exist. Some authors describe subharmonic response as an oscillation with 
exactly the frequency of the linearized problem excited by a force with three 
times this frequency. At least in the case of the Duffing equation, such an 
oscillation cannot occur. 

Interesting conclusions can also be drawn from relation (33), which de- 
termines the second Fourier coefficient of the subharmonic vibration as a 
function of A1;3. For A1;3=0 relation (28) is satisfied identically while relation 
(29) reduces to relation (8) for the “harmonic” case. Hence the subharmonic 
vibration results through bifurcation from the harmonic one. This takes place for 


A, =A — $4643. (36) 


In addition we find for A1/;3=0 
A 
dA;/d(w*) = — ry > 0 (37) 


and this determines the direction in which the new branch leaves the original 


one. 
Fig. 8 indicates the nature of the curves for A1/3 and | Aa| for both a 


hard and a soft spring. 
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Our calculations are accurate only for values of A1,3 that are small in com- 
parison with those for Aj, i.e., in the neighborhood of the bifurcation point. 
More interesting from the point of view of the applications are the solutions 
for which Ax; is large in comparison with A, i.e., those for which the sub- 
harmonic component dominates. It is very plausible that these solutions will, 
like their harmonic counterparts, lie in the neighborhood of those for the free 
oscillation if the force amplitude F is small. Of course, the response curve for 
the free oscillation in question would be the one beginning at w=3+/a and 
not that beginning at w= Va. 














\ 
Ay Ay 
a «= * = Ne 
IAI IA! 
Ma “ Wa ww 
Fic. 8 


If the subharmonic component dominates sufficiently so that the resultant 
solution x(@) is monotone in a half period from the maximum to the minimum, 
then it is possible to apply the method of Rauscher to obtain it. 

There are cases which are extreme in another sense, i.e., those in which 
the dominating harmonic component has a frequency that is a very small 
fraction (instead of 1/3) of the impressed frequency. In cases involving non- 
linear damping such phenomena occur—subharmonic oscillations up to the 
200th have been observed, apparently. 

9. Other periodic solutions. The harmonic and subharmonic oscillations 
so far considered do not exhaust the possibilities of periodic solutions. There 
exists a great variety of other periodic solutions which arise near w= V/a/n, 
n=2, 3, - +--+, in contrast with the subharmonic oscillations which arise near 
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w=n/a. Among them are periodic solutions in which the harmonic compo- 
nent of frequency mw is dominant. Since the fundamental frequency is w in 
such oscillations, they are not exact counterparts of the subharmonic ones; 
nevertheless they might be called superharmonic oscillations. They can be ob- 
tained by a perturbation method beginning with the solution of the linearized 
problem (@=0) x=A, cos wi+A, cos nwt for w= Va/n, A,=n?/(n?—1)F, 
while A, is prescribed arbitrarily. Such solutions have been investigated theo- 
retically, cf. [10, 12, 13], by a different method. They do not appear to have 
been observed experimentally. 
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ON THE DEFLECTION OF ANISOTROPIC 
THIN PLATES* 


BY 


VLADIMIR MORKOVIN 
Brown University 


1. Introduction. The elastic medium considered in this paper! is assumed 
to have at each point at least one plane of elastic symmetry parallel to the 
neutral plane of the plate, which is taken as the xy plane. The theory of bend- 
ing of thin plates possessing this type of anisotropy has been developed 
mainly by Boussinesq,? Voigt,* and Lechnitzky.‘ It is based on the usual as- 
sumptions of thin plate theory, leading to the following relations which are 
valid throughout the (small) thickness 2/ of the plate 





cae 0, (1.1) 
Ow Ow 
“=—Z ) v= —s—-» (1.2) 
Ox oy 


where w denotes the (small) deflection of the neutral surface of the plate, 
while u and v are the displacements in the x- and y-directions respectively. 

If the equilibrium conditions are used together with Generalized Hooke’s 
Law, (1.1), and (1.2), it is possible to express® the quantities characterizing 
the state of stress in terms of the partial derivatives of the deflection w. Then 
by considerations of equilibrium w is found to satisfy the differential equation: 











aw a4w he d4w 
bi3 — + 3bi6 - - + 2(di2 + des) =a a 3b26 
dx Ox*dy Ox*dy? dxdy* 
O4w 3q(x, ¥) 
pec, 0 
oy* 2h* 


where };; are constants depending upon the elastic material and g(x, y) is 
the normal load per unit area on the upper face, the lower face being free. 
The problem, therefore, is to find a solution w(x, y) of (1.3) which satisfies 


* Received Nov. 18, 1942. 

1 The author wishes to express his appreciation to Professor I. S. Sokolnikoff, who proposed 
the problem and suggested the manner of approach. 

2M. J. Boussinesq, Journal de Mathématiques (Liouville) Ser. 3, 5, 329 (1879). 

3 W. Voigt, Kompendium der Theoretischen Physik, Leipzig, 1895, pp. 436-455. 

4S. G. Lechnitzky, Applied Mathematics and Mechanics, Leningrad, 2, 181-210, (1938). 

5 An exposition of the theory is given in I. S. Sokolnikoff’s Mathematical Theory of Elas- 
ticity, (mimeographed, Brown University, 1941), pp. 319-329. His notation is adopted in the 


present discussion, 
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the prescribed conditions at the boundary (edge) Co of the neutral surface of 
the plate. This general problem has not been solved even in the case of very 
simple boundaries. Some success was achieved by various writers‘ investigat- 
ing special shapes of plates made of material with three mutually perpendicu- 
lar planes of elastic symmetry (orthotropic material). In the present paper a 
general method of solution is indicated when the boundary d, is an analytic 
curve. Detailed solution, illustrating the procedure, is carried out in the case 
of a clamped elliptic plate and a polynomial loading function g(x, y). The 
complicated form of the solutions obtained is due to the generality of the 
problem. For a given material the values of the parameters );; in (1.3) are 
specified and the computations become relatively simple. 

2. Preliminary considerations. Utilizing the fact that the potential energy 
of any realizable state of stress is always nonnegative, Lechnitzky’ proved 
that the roots yu, of the characteristic equation 


bir + 3dis w+ 2(di2 + bes) pu? + 3boe p* + bee uf = 0 (2.1) 
corresponding to (1.3), must be conjugate complex numbers, say, 
Mi = a + 183; Ms = A, = a — 18;; B, ~ 0; 


2.2 
2 = a2 — 182; Be ¥ 0. tnd 


:~) 


Me = a2 + 189; Ma = 
Consider the complex variables z; and 22, related to x and y as follows: 
Ze = et wey = (x + ary) + (Bey), k=1,2. (2.3) 


This formula specializes to zo, the complex variable in the original xy plane 
when uo =7. The relationship between the complex plane of zo and that of 2 
Or 22 1S 
Ze = pr20 + Qiko, k=1,2, (2.4) 
where 
Pc = 4(1 - ipx) and qt. = 3(1 _ ifi,). (2 .4a) 


Geometrically, (2.4) corresponds to affine transformations of the 2 plane. 
Since 6,0, these transformations possess inverses, given by: 


1 
Zo = Bs (PiZe — Queen), k = 1, 2. (2.5) 
k 


Let wo designate a particular integral of (1.3). Lechnitzky® has shown that 
the most general solution of (1.3) can be expressed in the form: 


* For instance, M. T. Huber, Probleme der Statik technisch wichtiger orthotroper Platten, 
Warszawa, 1929; W. McDaniels, thesis, University of Wisconsin, 1940; S. G. Lechnitzky, 
loc. cit., under assumption of vanishing normal load g(x, y). 

7S. G. Lechnitzky, loc. cit. 8S. G. Lechnitzky, loc. cit. 
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w = fi(zi) + falBs) + fo(ze) + fo(Z2) + wo, (2.6) 


where f;,(z,) are arbitrary functions of z;, and f(z) are their conjugates. These 
functions are analytic in the regions which correspond through (2.4) to the 
region occupied by the plate. If the roots 4; and we in (2.2) are not distinct, 
the solution w assumes the form 


w = 2ifi(2r) + 2rfi(21) + giles) + gilZs) + wo. (2.7) 


When the region in the xy plane is simply-connected, these analytic func- 
tions are also single-valued. The present discussion will be confined to simply- 
connected regions bounded by analytic curves Co. Let the parametric repre- 
sentation of the boundary be: 


Co: z= x+ ty = F,(2) + iF,(2), 035% < Ze: (2.8) 


The period of the functions F,(¢) is 27; their Fourier coefficients are 
1 2r 1 2 f 
day = —f F,(t) cos nt dt; ens = —{ F,(t) sin nt dt. (2.8a) 
WT 0 us 0 


The boundaries of the corresponding regions in 2, planes, by virtue of (2.3), 


are 
Ci: + F,(t) oo ul 2(t), k = iF Ye (2.9) 


It should be noted that the homogeneous deformations (2.4) and hence 
the contours C; are determined at the outset by the anisotropy of the plate. 
The controlling parameters are five in number, namely by, die, (b12+ Dee), 
bes, bee, or bi1, a1, Bi, G2, Bo. 

3. Outline of procedure. The form of the solution (2.6) suggests that some 
simplification could be achieved by considering the boundary conditions di- 
rectly in the z, planes. The boundary conditions usually consist of two func- 
tional relations between the deflection w and its partial derivatives of at most 
third order on Cp 


ow Ow 0*w 
G, | w,—» —>-+-+,—]=0, *s = {, 2. (3.1) 
Ox dy dy® 
If one regards z, and 2, as the independent variables, then 
fe) 0 + fe) re) ts) te) (3.2) 
a ee oe teed . ie-eadie dl. =a J. 
ax Os, a% ay as, | ae 


Hence, the boundary conditions (3.1) can be formulated along either C; or C2: 


ik dw dw aw 
G,| w,—» —»+++,—])=0, ng, k = 1, 2: (3.3) 


OZ; OZ, 073 
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It is desirable to express the boundary conditions in terms of a single 
variable. To effect this, suitably chosen neighborhoods of contours C; in 2, 
planes are mapped conformally into some neighborhoods of circumferences <x 
of unit radius in new complex ¢; planes in such a way that points on C;, cor- 
respond to points on yx. Let g, represent the polar angle in the ¢; planes. 
Identical values of @ and ¢g2, 0S¢2<27, can be made to correspond to one 
and the same value of the parameter ¢ on Co, 0S¢<2z, by the proper choice 
of the mapping functions 


z= wel(€x), k= i, 2. (3.4) 
Laurent expansions for w,({%) are found in terms of the Fourier coefficients 


of F,(t). Let the common value of ¢; and {2 on yx be denoted by o=e*; then 
by virtue of the relations 





0 1 re) 0 1 rs) 
en ae oe ee (3.5) 
Oz, wel(Se) Oe OZ, we(Fx) Oe 
the boundary conditions (3.2).can be expressed in terms of ¢ alone: 
Ow ow Ow 
r,( w,—» —>-++,—)=0, #2 = 1, 2. (3.6) 
00 «0g 0c 


The deflection w may have the forms 


1(f1) + 1(F1) + b2(t2) + b2(F2) + wo, (3.7) 
w= o1(Fi)or($1) + w1($1)b1(F1) + oilSr) + ilFi) + wo, (3.8) 


WwW 


where 


O(n) = ful o(Sx)} and yultx) =ge { won(tx) } 


are undetermined functions, analytic and single-valued in some neighbor- 
hoods of y;. Accordingly one has 


+00 +00 
alte) = Do Ameves Velde) = DO Bude - (3.9) 


Expanding I’, in (3.6) into powers of ¢ and equating coefficients of like powers 
to zero, one may expect to obtain recursion formulae for the determination 
of the coefficients A», and B,,, in (3.9). This, in essence, solves the problem. 

4. The mapping functions. The mapping functions o,({;), if they exist, 
must be analytic and one-valued at least in narrow rings around the circum- 
ferences yx where they can be represented by convergent Laurent series, say 


+00 
t= >> Dintt, k=1,2 (4.1) 


m=—w 


® This will be recognized as an extension of the scheme of N. I. Mushelisvili. See Mathe- 
matische Annalen, 107, 282-312, (1932). 
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The correspondence between points on the curves C;, and yy; is expressed 
through the equation 


+0 +2 
Cy: 2. >= y Dimox = } Dim(cos Mer + i sin m,). (4.2) 
The parametric representations of C; are: 
Cut te = (dor + uedor) + D> { (dni + wedne) cos mt + (€n1 + pxena) sin nt}. (4.3) 
n=1 


Letting identical values of g: and ¢g2 correspond to the same value of the 
parameter #, yg, =#, one can solve for the coefficients Di» from (4.2) and (4.3): 


Din = 4{ (dni + apdne + Benz) — i(€ni + axenz — Bxdn2) i, 


; (4.4) 
Dy _ AY (dni + apdne — Bren) + i(€n1 + Aren2 + Bidno) } ’ n a 0. 
It will be shown presently that 
Tim | Din] < 1, lim ~/| Din] < 1, (4.5) 


n-~0e no 


so that series (4.1) with coefficients (4.4) actually converge in some rings 
around y; and the foregoing formal steps are justified. The proof rests on the 


following theorem:?° 
The necessary and sufficient condition that the periodic function f(z) defined 


by the series 


c] 
tao + >> (a, cos nz + b, sin nz) 


n=l 


be analytic in some strip parallel to and containing the real axis 1s 


lim */(a2 + 62) <1 or lim | a, + ib,] <1. (4.6) 


ne n—~2o 





The functions F,(t) = 3do.+>_7.;(dns cos nt+e,, sin nt), s=1, 2, are by hy- 
pothesis (real) analytic functions. Hence the functions defined by the series 


1 tC ) 
> DX { (dar + ardne + Brenz) cos nt — (nr + oxen: — Budn2) sin nt}, 
in (4.7) 
2 z. { (dai + apdne — Beene) cos nt + (€n1 + Aren2 + Bidno) sin nt} ’ 
n=1 


are also (real) analytic. Replacing / by the complex variable 2 which reduces 
to ¢ on the real axis, one is led to complex functions, satisfying the conditions 





10 This is a slight modification of a theorem found on pp. 125-127 of De La Vallée Poussin’s, 
On the approximation of functions of a real variable and on quast-analytic functions, Rice Institute 
pamphlet, 1925. 


1943] ON THE DEFLECTION OF ANISOTROPIC THIN PLATES 121 


of the preceding theorem. But since (4.6) applied to the functions (4.7) im- 
plies (4.5), the existence of regions of convergence of (4.1) containing ; is 
shown. 

These regions can be narrowed to rings about" yy, in which the deriva- 
tives of the mapping functions have no zeros so that the functions themselves 
are simple (schlicht). The representations of the narrowed regions on the cor- 
responding bands about C, are therefore one to one and conformal as desired. 
The foregoing discussion can be extended to any number of planes 2,. If, for 
instance, one sets “9 =17, one obtains a function wo({») mapping conformally 
the original curve Cp together with a surrounding band into yo, the circum- 
ference of unit radius in a new fo plane, and its neighborhood. 

5. Mapping of elliptic regions. The procedure described in the preceding 
section constitutes a practical scheme for the determination of the mapping 
functions w,;(¢;) once the shape of the plate is known. In the case of the ellip- 


tic boundary 
Co: 2% = acost-+ ib sin t, (5.1) 


formulae (4.1) and (4.4) yield the mapping function 


Ri 
Ze = or(fe) = (te +—), (5.2) 
ve 
where 
Tk 
R, = —; re = 3(a + ip,d); Sze = 3(a@ — ip,b). 
Sk 


It follows from elementary considerations that the deformed contours in 2, 
planes 
Ce: 2 = acosg+ u,b sing (5.3) 


are ellipses whose major axes make angles 6; with the real x, axes, where 


9 2 
tan (26,) = zone? = am { 4ryse } ° (5 4) 
a’ + b*(ag — Bj) 

It also follows from (5.2) that the families of circles concentric with y; in the 
¢, planes correspond to the families of confocal ellipses with foci at 2\/ns, 
in the z; planes. In particular, to the circumferences y; correspond the ellip- 
tical contours (5.3) and to the smaller? circumferences with radii V|Rzl, the 
degenerate ellipses of the families, i.e., the double segments joining the foci. 

These degenerate ellipses can be represented in the form 








11 », may well be the outer boundary of these annuli when the derivatives F! (¢), s=1, 2, 
vanish simultaneously. Thus the class of curves to which this discussion is applicable is some- 
what larger than the class of analytic functions as usually defined. 

12 If uw; and yw. are defined to be the roots whose imaginary parts are positive, 8, >0, then 
V/| R,| <1. A 
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ze = 2\/(risy) cos (yp — 5x) (5.5) 


where 6; are the amplitudes of ~/R;, the ¢; plane images of the foci 2\/7s,. 
Since the derivatives w/ (¢,) vanish only at \/ Rx, one may choose for the rings 
in which the mapping functions are to be schlicht the annuli 15$;%<~/|R,|. 
The corresponding regions in z; planes comprise the interiors of the ellipses C; 
made doubly-connected by the introduction of slits (5.5) between the two foci. 

From (5.5) it is clear that every point 2; on the slit corresponds to two 
distinct points, | R,|e” and VR, | ei, on the inner boundary of the 
above ring. Since the analytic functions f;(z.) are single valued inside C;, 
the transformed functions ¢;(¢,;) must assume identical values at the two 
points. Hence, 


+o +00 po ae 
Zz. AmkN | Ri \™ eme = y A mk | Ri. im (2be—e) | (5 .6) 


m=—c m=—oo 


so that 
A_,.k = RiAn, n> 0, 


and finally 
R n 
k ro 
). (5.3) 
oe 
The question arises whether the restriction (5.6) is sufficient to insure that 
the functions F;(z,) obtained by the inverse transformation 


1 ——— a 
a= {2% + \ (22 — 4rxsx) } (5.8) 
2s% 








C) , a R. 
); vi(S) = > Bak («i + - 


n=0 


otd= > Aua( ti + 
n=0 


from the functions ¢;(¢;) in (5.7) are single-valued in the full, simply-con- 
nected interiors of the ellipses C,. An affirmative answer is reached after ap- 
plication of Schwarz’s Reflection Principle to the potentially different branches 
of F,(z,) on the opposite sides of the slits. Therefore, if the unknown coeff- 
cients in (5.7) are determined from similarly transformed boundary condi- 
tions, the resulting functions will correspond to the solutions of the original 
problem. 
6. Uniqueness of solutions. If in (2.7) the functions f;(2:) and gi(z:) were 
replaced by fi(z:) +f#*(z:) and g:(2:) +g:*(21) where 
Re{z:f# (21) + g#(z:)} = 0, (6.1) 
the value of w would remain unaltered. The extent of arbitrariness in the 
choice of fi(z1) and gi(z:) is determined next by finding the most general 
analytic functions f*(21) =ui+im, and gi*(z:) =u2+iv2 satisfying (6.1), or the 


equivalent condition 
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Xu, + V1 tw =0. * (6.2) 


Utilizing the fact that m1, v1, ue, and ve are harmonic functions, conjugate in 
pairs, one readily finds from (6.2) that 


Ou, Ov) 
and therefore 
fi (21) = Arizi + (Ae + AAs), (6 .3a) 


where A; are arbitrary real constants. Substituting this expression into (6.1) 
one arrives at 
gir (21) = (Ae = id3)21 4 1d4. (6 .3b) 
The foregoing functions are transformed by (5.2) into 


Ry 


*(¢1) 


miss(t + )+a4 13; 


$1 


(6.4) 


R, ; 
— (Ao — tA3)51 («: + ~) + 1r4. 


1 


vir ($1) 


It is desirable to eliminate the arbitrariness inherent in the form of solutions 
by choosing the parameters \; so as to simplify the expansions (5.7) of $1(f1) 
and y,(¢;). If 


—2 Im {Ani} , : 
A1 = a do + 1X3 = — 2A 015 Ng =--—2 Im {Bo}; 
a + bB, 


- . R. " R n 

¥ 4n(# +—) and ¥ Bul 1 +—), (6.5) 
n=0 $f 

with Ay, =A, and Bo, = Bo. That the solutions of the form (6.5) are uniquely 


determined by the boundary conditions is shown in Section 8. 
A similar argument shows that one needs to consider only the functions 


0 = R. rc ) J R: 
oi(f1) = D An(# + =); o2(f2) = ¥ An +=), (6.6) 
$y 


n=0 n=2 2 


with Ao, =Aoi and Ax =A, when the roots 4; and pe are distinct. 
7. Boundary conditions. Let @ designate the angle between the outward 
normal to the boundary Cp» and x axis. Then, 
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1 
sin @ = — — (dzo + do) =—— (ift,dz, —_ i,dZ,), 
2ds 2dsBx . 
(7.1) 
1 1 
cos § = (dzp — d%) = — (dz, — dx), 
2ds 2dsB 


where dzp and ds are differentials along Cy and dz, the corresponding (depend- 
ent) differentials in the z; planes. The equations (7.1) and (3.2), together with 
the mapping functions (3.4) make it possible to formulate the boundary con- 
ditions in terms of ¢ on the circumferences y;. In the case of a clamped plate 


one has the conditions: 


Ow fs) ; re] : 
w(x, y) = 0 — =(cos 0 — + sin 0 — } w(x, y) = 0, (7.2) 


W(X . 2 
On Ox oy 


on Co, which are equivalent to the vanishing, on yx, of the expressions* 


wio,o) = 0, (7 .3a) 
1 = { 5 ri 1 a; (a) , } re) 
-_-— — o( Didi + OG) — ——— 29? 
Bio dst \ 4 td a wi (c) Pid Oo 
J we (c) 1 \9 " “: 
+ 4¢ a ce, 2D igh + — (pibi + gigk) ae w(e, g) = 0, (7 .3b) 
{ w; (oc) o f ag 


where do is dependent on dzp. Obviously, to the part of w(¢, ¢) obtained by 
transformation of f;(z;) and gi(z:) one applies the operator in (7.3b) with k =2. 
It is expedient to make use of the mapping function zo =wo({o) for the trans- 
formation of the particular integral wo. When k=0 the operator in (7.3b) 


reduces to 








it do 0 i 2 ; a 
—-— —[g — + — =) wi, Gc), (7 .4) 
a ds Oo a Oa 
and (5.2) has the form 
. 1 a+b ‘ 
Zo = wWo(Co) = 3(a = b) ( RF + —), ke = ’ (7.5) 
fo o= b 


If the normal load g(x, y) in (1.3) is expressed as a polynomial" in x and y 
of degree (m:—4), wo is a polynomial of degree m. For instance, when the load 
is linear 


2h® 

q(x, y) = — (co + 1x + Coy), (7.6) 
3 

where c; are real parameters, one can take for wo, 

13 Tt is noted that, as a consequence of (3.2) and (3.4), o and @ are treated as independent 


variables in the partial differentiation process, whereas the result can be simplified by o=1/o, 


an identity on y,. 
4 This may be considered as an approximation to more complicated loading functions. 
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Co 
Wo(x, y) = —— x! x5 + v5, ‘| 
, i 4! by, 5S! by, 5! Boo ; ( ) 








The following abbreviations referring to values on the original curve Co but 


expressed on ,; are introduced: 
ny 1 
Ww=- K(c) = = ¥ (no + bn =). 
0 o”™ 
ds OW 


ig — — 
-dao On 


(7.8) 


II 


ny , 1 
“ws L(o) — ¥ (ono + am -). 
0 


o a 


One readily finds the coefficients a, and 6, for the special case of the linear 


load applied to an elliptic plate: 


4ba*co 5ab (= _ €2b° 
a So. el a5 = ice ee ae 
244! b1, 25 5! bis boo 
ag = 3a4, ag = 4a; a3 = 5ds, a, = 10az5, (7.9) 
a a . 
bom = — Gem} (boma41 + bom+1) a So (Gom+41 > Gom+1)- 
4b 5b 


The preceding transformation of the particular integral wo is valid whether 
the deflection w has the form (2.6) or (2.7). However, the parts of w involving 
the unknown functions f;(z,) and g;(z,) have to be treated somewhat differ- 
ently in the cases of equal and distinct roots. To avoid repetition only the 
first case is discussed here. 

8. The case of equal roots. Substituting'® (3.8) into the boundary condi- 
tions (7.3a) and (7.3b) and simplifying gives: 


1 1 
a( oto) + o(o)6(—) + ¥(o) + 1(—) = K(o), (8. 1a) 


1 1 
[ (. -) ¢'(c) + w'(a)g (-) + vo) | 


w'(c) 1 
+—[o 250 + = (08 + | 
o 


1 
[oor (2)+2(Doote(L)]=20. ew 


16 In this section the subscripts are omitted since the distinction between 2; and 22 planes 


is not involved in the discussion. 
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In order to avoid a system of an infinite number of equations in infinitely 
many unknowns the following device is used. The condition (8.1b) is replaced 
by a combination of (8.1b) and (8.1a). The first boundary condition (8.1a) 
is an identity in o and will therefore yield an equality upon differentiation 
with respect to o. This equation is multiplied by 


g ( ow'(c) : 
| laa 


and added to (8.1b). In the resulting identity, one solves for the expression 
involving the unknown functions and obtains* 





1 1 
s(-) a¢’(c) + ow'(a)d (—) + o’(c) 
o 


” s{-5 — De 4 —| 
2 1 i] 1 
Lae (—) + gow'(c) pow'(c) + hk. a’ (-) 


o 





o o 


1 


ot 


Lic) {1\ 6K") "\e : 
{=e (=) + = (eu'(or2h + —~— (08 + a0) (8.2) 


o oc 


When the plate is elliptic, (8.1a) and (8.2) read: 


1 os oo R™ R ra) a 1 a: 
§ (— + Re) ; ’ A m (om + ) + 5 (« + -) ; A m ( + Rvo™) 
Co , o”™ o i o™ 


oo R™ w# _ 1 _ 
+> Bn (om + — + > B,, (— + Rvo™) 
o”™ 0 


0 





m1 1 
=> (520 + Bn -), (8.3) 
0 o”™ 
1 _ : z* R ae ] = 
s(— i Re) > mAm (om - —) aa s(o — —) A m (= + Rea") 
o 1 _ o 1 ” i 
cs) R™ 1 po? g | 
Pees a 
" X " (« =) oB(a + 3B) F — o*/k? 1 — 1/k?o? 
1 - m1 Ian 
{s(= _ Re) be (ano + =) + 4} (d — ifia)e 
o 0 o™ 


1 mt 1 
+—(b+ ijza) } >» m (dno — bn =)]. (8 .4a) 
. o 0 o 


%6 The expressions in the denominators, 1/8(pow'(c)+9/o @'(1/c)) and 1/8(p/o @'(1/c) 
+gow'(c)), never vanish, since pp —qq=8 £0. 
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For a circular plate a=b, 1/k® vanishes, and the last equation is already ex- 
panded in series of powers of a. Otherwise, the expansion of the right hand 
member is infinite, since 


1) fh pee SE ws 


o \ 1 — s?/k? 1 — 1/k?o? 1-0 \kR o 1-0 \ko 


Combining the terms 3(1/a— Ro) and 34((b—ifa)o+1/o(b+iga)) with the 
expansion (8.5), one simplifies the right-hand member of (8.4a): 


ny 1 se 1 
| (o-+und) > (ano"-+a0 =) i(jia— pb) >» m (sn0"—5 ~)| 
0 g” 0 o”™ 


28 a+b) 


2ab | (Ef denotes 5) 1) (=) 
B(a?— 6) Pq - am —M0y)o amTM0m = k 


+ pa( | (an-t-mbnom-+ (an — mb) -\) >(=) |. (8 .4b) 
o” o 


0 1 


Let Q, denote the coefficient of o” in (8.4b). This coefficient is evaluated by 
means of relations similar to (7.9). It is found that Q, vanishes for n>. 
Recalling that Ap =0, A;= Aj, and Byo= Bo, one concludes from the rela- 
tions between the coefficients of o° in (8.4a) and (8.3): 
Qo 1 


11+P 
A, = cece Bo os — bo — 
(1— P)(s+5) 2 41=-P 





Qo, where P= RR. (8.6) 


The system of equations between the coefficients of o*, »>1, 


5Ansi + SRR"*"Any: + B, + RB, = 5, — SRAg-1 — SR*™'An-1, 
(8.7) 
3(m + 1)Angs — SRR**'A,41 + mB, = 0, — SR(m — 1)An_y — SR™A-1, 


and of the conjugate equations has a non-vanishing determinant: 
s$Dn4i(P) =s3{1—(n+1)?P*+2n(n+2) P*+!—(n+1)?P+24+ Pet?) m>1, 


because P <1. The system, therefore, possesses a unique solution given by 


the expressions: 


1 —_ —_ 
Anyi = a [On(1 — (m + 1)P™ + nP™+1) —0,R(n — (n + 1)P + P**) 
SL/n4+1 
— nb,(1 — Pt!) + n(n + 1)5,R™(1 — P) + SRD,(P)Az-1 


+ sR""E,(P)A,-1], (8.8) 
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—1 
B, = om [O,{1 — (wn + 2)P"t! + (m + 1)P**} 


—O,R"{ (nm + 1) — (n+ 2)P + P»t*} — (nm + 1)b,(1 — Pt?) 
+ b,R”} (n +1)? — n(n + 2)P = prt?} 
_ $RGni1(P)An-1 + sR" Eni 1(P)An-_1], (8.9) 


where the following abreviations are used: 


E,(P)= | (n—1) —2nP+(n+1)P?+(n+1)P*—2nP™*!+ (n—1)P**?}, 
Gnii(P) = | 2—n(n+1)P*-!+ (n—1)(n+2)P*+(n—1)(n+2)P™*! 
—n(n+1)P™t?42P2nt1} | 


One can calculate any number of coefficients in the expansions of $(¢) and 
¥(¢) from those already known, Ao» and Aj, and the recurrence formulae (8.8) 
and (8.9). If (8.8) is rewritten for n>, in the form 


Anti 1 _— 7 An-1 
ae = RDA?) + —R+1E,(P) =I (8.10) 
ye Dnii(P) § 


4in—l 


One sees that limy +0| Any1/A n-1| =|R|. Consequently, ¢(¢) and ¥(¢) con- 
verge and represent analytic functions in the ring | R|v2< |e] <|R| /2, 
which extends beyond the ring corresponding to the interior of the ellipse Co. 
The expansions are finite if A,-; and A, vanish for any m>™m,. In the case of 
the linear load one finds As =Asg=0, Bs¥0, so that the solution is a poly- 
nomial of fifth degree. Specific examples indicate that the solution is a poly- 
nomial when the loading function g(x, y) is a polynomial. 

9. Conclusion. The discussion in the case of distinct roots is similar to 
that of Section 8. For details and results the reader is referred to the author’s 
doctoral thesis.” 

When si=42=7, the solution reduces to the results for an isotropic 
clamped elliptic plate bent by uniformly varying load.'* Another impor- 
tant special case is that of an orthotropic elliptic plate bent by a linear load, 
for which the solution is new. There are two possibilities obtained by setting 
either 4; =781, M2 =782, or wy =Q,+781, Ue = —a,+78;. The calculations in these 
and other specific cases are comparatively simple. It is clear that by specify- 
ing the values of some or all of the parameters characterizing the material 


17 On the deflection of anisotropic thin plates, University of Wisconsin, 1942. 
18 See A. E. H. Love, Treatise on the Mathematical Theory of Elasticity, Cambridge, 1927, 


p. 486. 
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of the plate, its shape, and the type of loading, a considerable simplification 
is achieved. 

The foregoing method can be extended to the cases of (a) different bound- 
ary conditions, (b) different shapes of plates. Many of the devices and formu- 
lae, such as those of Sections 4, 6, and 7, remain valid. The extension to the 
interesting case of infinite doubly-connected regions involves little more than 
the additional feature of determining the character of the multiple-valued 
functions f;(z,) from the fact that the deflection w remains single-valued and 
continuous. 
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SYMMETRICAL JOUKOWSKY AIRFOILS 
IN SHEAR FLOW{ 


BY 
HSUE-SHEN TSIEN 
California Institute of Technology 


1. Problem. The usual two-dimensional theory of airfoils assumes a uni- 
form velocity for points far from the airfoil. There are many applications 
where this condition is not satisfied. For example near the ground there is a 
large vertical velocity gradient, and therefore a first approximation to the 
problem can be obtained by assuming a linear velocity distribution. This has, 
in fact, been done by H. v. Sanden! in connection with O. Lilienthal’s experi- 
ments in natural wind. However, v. Sanden used a numerical method of in- 
tegrating the differential equation and carried out the calculation only for a 
wedge-shaped body. Th. von Karmén* suggested to the author to take up 
this problem again in order to develop a more complete theory. Hence in the 


y 





Fic. 1. Body in a Shear Flow. 


t Received Dec. 27, 1942. 
1y. Sanden, H., Uber den Auftrieb im natiirlichen Winde, Zeitschrift f. Math. u. Phys., 


61, 225, (1912). 
* The author wishes to thank Dr. von K4rm4n for suggesting the problem and his kind 


interest during the course of the work. 
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first part of this paper, a generalization of the well-known Blasius theorem 
for calculating aerodynamic forces acting on an airfoil is given. Then the 
result is applied to the case of symmetrical Joukowsky airfoils and the final 
data are given in a number of tables and graphs. 

2. Method of solution. By setting up the problem as shown in Fig. 1, the 
velocity distribution far from the airfoil is given' by 


v(t + x2) (1) 


c 


u 


9=0 


where Up» is the undisturbed flow velocity along the x-axis and K is the non- 
dimensional velocity gradient of the undisturbed flow. c is some dimension 
of the body immersed in the stream, e.g., the chord of the airfoil. Then the 
vorticity at locations far from the airfoil can be calculated from Eq. (1) as 
ov Ou UoK/ (2) 
——_«as «= SE a= 0 Cc 
Ox ody 
which is a constant. However, in the flow of non-viscous incompressible fluid, 
the vorticity is associated with the fluid and maintains its strength. Consider 
the field of flow starting from the far left, where the vorticity is constant and 
equal to — U,K/c. This value of vorticity is carried with the fluid over the 
whole field of flow. Therefore, the flow problem on hand is one with constant 
vorticity distribution. 
To satisfy the equation of continuity, 


Ou av 
—+—=9 (3) 
Ox dy 
the stream function y is introduced. It is defined by 
dy oy 


vs re) 7 ~ Ox (4) 


u 


u, v are the components of velocity in the x-direction and the y-direction. 
Due to constant vorticity distribution, the vorticity equation is 
Ov «Ou 


—=——«=«a — U,K/c. 5 
ax dy oK/ (5) 


By using Eq. (4), Eq. (5) can be written as 


re] 3? 
a UoK/c. (6) 
dy? 


Ox? 


Therefore the flow problem is reduced to that of solving Eq. (6). 
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The stream function Wo of the undisturbed flow given by Eq. (1) is 


K ¥? 
Yo = U0(y + 7 =) (7) 


c 


which can be easily verified by means of Eq. (4). The mathematical problem 
can be considerably simplified by introducing the stream function ¥; due to 
the presence of the body defined as 





Y=vWty. (8) 
By substituting this expression for y into Eq. (6), the equation for y; is simply 
0? 0 
i <é (9) 
dx? dy? 


This is the Laplace equation. Therefore any solution of the Laplace equation 
x Q5 10 45 


-/0 “GS 0 











TO ye? */ 






Y~ ay rye t( tax -2) 

















Fic. 2. Source in a Shear Flow. 


combined with Yo will satisfy Eq. (6). For example, we can combine Wo with 


a source 
¥1 = Uobd! (10) 


or a vortex 
¥i = Uoad logr (11) 


as shown in Figs. 2, 3, and 4. Here 


V —— 
6 = tan-!'—,; y= VJ/x? + y’, (12) 
x 





1943] SYMMETRICAL JOUKOWSKY AIRFOILS IN SHEAR FLOW 133 


It is interesting to notice that in the case of a source the zero stream line, 
which forms the walls of the “half body” when K =0, is no longer symmetrical 
with respect to the flow direction. The velocity, and hence the pressure, along 


os 10 IS 











ee Y= 2y + y*+ 48 (ton ¥-n) 
-O¥ 


ins --05 


-06—_ 


oe S) + -10 
-07- 
TO y-- (3 + 1) 


Fic. 3. Source in a Shear Flow. 














Fic. 4. Vortex in a Shear Flow. 
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these two branches are also different. If the flow field within the zero stream 
line is replaced by a solid body, such lack of balance may cause the resultant 
pressure force to become infinitely large with the infinitely long solid bound- 
ary. This surmise will be verified later when the resultant force and moment 
are calculated. 

If the boundary of the solid body is given, the flow problem indicated in 
Fig. 1 subjects y to following conditions: 

(i) the disturbance velocity due to y must vanish at points far from the 


body so that Eq. (1) is satisfied; 
(ii) the normal component of the disturbance velocity at the surface of 


the body must equal the negative of that due to Wp so that the resultant nor- 


mal velocity is equal to zero. 
3. Shear flow over a circular cylinder. To illustrate this method before 


solving the more complicated case of an airfoil, the flow over a circular cyl- 
inder will be investigated first. If the center of the circle is located at the 
origin, and the radius of the circle is c/2, the undisturbed velocity U due to 


Wo at the surface of the cylinder is 


K 
U= (1+ = sin s) (13) 


where @ is given by Eq. (12). The normal component U, of the velocity U is 


oa 


U, = U cos 0 = Ug (cos 6+ = sin 28). (14) 


Therefore, the normal component of the disturbance velocity due to ¥; at the 
surface of the circular cylinder must be equal to — U,, or 


1 oy; K 
(— =) = — Uy (cos 6+ — sin 28) (15) 
r 060 /rmc/2 4 


where 7 is given by Eq. (12). 
On the other hand, the solution of Eq. (9) that will give vanishing dis- 
turbance velocities at points far from the origin is 
* 1 
¥i = Uo E log r + bo9 + b> (a, cos nO + b, sin n@) =|. (16) 
n=1 ia 
The a, and }, are undetermined coefficients. By substituting Eq. (16) into 
Eq. (15), one has immediately 


n=-(5) 


a= (=) =. (17) 
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All the other coefficients vanish. Therefore, the resultant stream function is 





= - v,| ( “") sin 0 + (Saints “cos 26) |. (18 
¥Y=vythn= [(> <) sin +> ~ sin +o )} (18) 


At the surface of the cylinder, r=c/2, Eq. (18) reduces to 


(19) 








Fic, 5. Circular Cylinder in a Shear Flow with K =2. 


which is a constant, verifying the boundary condition of the problem. The 
stagnation point on the cylinder can be calculated by means of Eq. (18). The 
condition is (@~/0r),..=0. This is satisfied at the point where 


1 OK? 
sing = =| 4/1+— +1]. (20) 
K 4 


Therefore, if —8/3<K <8/3, there will be two stagnation points on the cyl- 
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inder located on that half of the surface where the velocity is higher. If K 
has a value beyond this range, there will be four stagnation points on the 
surface. 

Fig. 5 shows the stream lines for the case of K =2. It is seen that there is 
an additional stagnation point in the flow field on the negative y-axis. 

4. Force and moment. In this section, the force and the moment acting 
on a body whose disturbance stream function y is given by Eq. (16) will be 
calculated. 

The pressure p in the fluid is related to the velocity components u, v and 
the density p by Euler’s equations 


Ou Ou Op 

pu Peg pv [ =— = = 
Ov dv Op 

aa af allies i 


By means of Eqs. (4), (7) and (16), the quantities on the left hand side of 
Eq. (21) can be calculated and then p can be obtained by integration. How- 
ever in the final calculation of force and moment, the pressure will be in- 
tegrated along a contour far away from the body; only terms up to 1/r? in p 
need to be considered. Furthermore, the calculation can be simplified to a 
certain extent by differentiating the first of Eq. (21) with respect to x and 
the second with respect to y and adding the resultant. Then by using Eggs. 


(3) and (5), 
Su’ ES  . (+ UK ~'}. (22) 


p Ox? dy? axdy cc |6Of? 





By substituting Eq. (16) into Eq. (22), the differential equation for p is then 


a 1f1 0 / op 1 0p 
a V’?p = — - —— ro) + —— amen 
p aL? -orX\ OF r? 00? 


2 1 2 2 K (1 
2U 5 | - = (ao + bo) + —{— (— ap cos 20 + bo sin 26) 
r c tr? 


I 


1 
a — (2a; cos 30 + 26; sin 36) 
fi 
1 
ob ; (6a, cos 40 + 6b2 sin 40) +--- \]. (23) 
7 


The appropriate solution for this non-homogeneous equation is evidently 
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— 1 2 2 K 1 . 
— p = 2Uo| — — (ao t+ bo) + — 4 — — (— ao cos 26 + do sin 26) 
4r? c 4 
1 
— ma (2a; cos 30 + 2b; sin 36) 
r 


1 
— —— (6a, cos 40 + 6b2 sin so) | 


12r? 


1 
+ A+ Aologr + Bod + — (Ai cos 6 + B; sin 6) 
Tr 


1 
+ — (Az cos 20 + Bz sin 20) +--- (24) 
72 


where the A’s and B’s are undetermined constants. Either the first or the sec- 
ond equation (21) can be used to determine Ao, Ai, -- + and Bo, By, Be, - - - 
The final result for the pressure p can be written as 


p of K K K , 
—=A+ Us| — a log r + — bo — — (— ao cos 260 + dg sin 26) 
p c c 2c 


x 
— — (a, cos 36 + b; sin 36) 
2cr 


4 


1 3 Ka, 3 Kb, , 
+ —{(= — bs) cos @ + | — ~~ as) sin o\ 
r 2 Cc 2 Cc 


1 ao + by K : 
+54- —_——— — — (az cos 40 + bz sin 46) 


r? 2 Cc 


Kaz Kbz 
4. (2 — — bs) cos 20 + 2(— + a) sin 26k id ‘| (25) 
Cc Cc 


Now by considering the pressure force and momentum of the fluid, the 
following relations can be obtained between the forces X, Y acting on the 
body and their moment M, about the origin? (Fig. 1). 


X= - f pay — f ou(uay — vdx) 
c c 


Y= f pdx =f po(udy — vdx) (26) 
Cc Cc 





My, = f p(xdx + ydy) — f p(vx — uy)(udy — vdx). 
c c 
2 See for example Glauert, H., Aerofoil and airscrew theory, Cambridge University Press, 
London, 1930, p. 80. 
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The integrals are taken along any closed curve C enclosing the body. If a 
circle with radius R is taken as the contour, then in the integrals of Eq. (26) 


ll 


x R cos 6 


y = R sin 0. 


7 


(27) 


Therefore, by using Eq. (25), 


al 2 o. 
~ f pdy == Rf p cos 6d = apU E wa a a| . (28) 
c 0 phe 


Furthermore, 


~ 


K 1 
udy — vdx = Uy E R? sin 206 + Rceosé@+ bg + ry (— a, sin 6 + 6b, cos @) 
" 


1 
+ R (— 2a2 sin 26 + 2b. cos 20) +--- |« (29) 


and 


K 1 
“= Us| — R sin 6 + 1 +— (ay sin 0 + by cos 0) 
c 


1 
+= (= ax sin 20-+ 8, cos 20) +--+ | (30) 


By combining Eqs. (29) and (30), the second term in the equation for the 
horizontal force X can be calculated as 


; 3K 
— f ou(udy —vdx) = — roe] 3b. - =a], (31) 
Cc 2C 
Finally the horizontal force is expressed as 
= 2 
X = — 2xpU obo. (32) 


Similarly, the vertical force Y or lift and the moment M, about the origin are 
obtained in the following forms: 


~- 
‘ 
Il 


- 
2apvi| a ih ~< 6) | (33) 
Cc 


2 K 
M, = 2rpU | onb — @ + a (3Rbo = | . (34) 
c 


Eqs. (32), (33) and (34) show that the force and moment on the body can be 
calculated in terms of the strength of the vortex, the source, the doublets 
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and the quadruplets in the disturbance stream function y. These equations 
can be regarded as the generalization of the well-known Blasius formulas. 
They reduce to the latter formulas if K =0. However, it should be noticed 
that if there is a source, i.e., b) #0, both the lift Y and the moment M, 
grow to infinite magnitude with R- . This confirms the surmise stated pre- 
viously. If the boundary of the body is closed, by) must vanish, and there can 
be no horizontal force or drag and the lift and moment will remain finite. 
v. Sanden! obtained a small drag force for his wedge shaped body. Evidently 
this is due to the unavoidable inaccuracies in his numerical method. 

5. Symmetrical Joukowsky airfoils. For the flow over a symmetrical 
Joukowsky airfoil, it is difficult to determine the disturbance stream function 
¥ directly. But, as seen from Eq. (9), Yi satisfies Laplace’s equation, and 
therefore conformal transformation can be used. It should be noticed, how- 
ever, that Wo does not satisfy Laplace’s equation; and therefore it does not 
allow the use of this transformation in the ordinary sense. 


id 
a ae 


i] 


C| 
a 2 


(a) 








>< 


z- PLANE 





a _ 
R : 
KK 


c 


ny 


¥ 
*%, 
Crags 





(6) 


Fic. 6. Joukowsky Transformation of a Circle into an Airfoil. 


Consider a circle of radius a in the ¢ plane (Fig. 6a). The transformation 


( ++) (35) 
ence ake f—e ate 





will transform the circle into a symmetrical Joukowsky airfoil in the z plane 
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(Fig. 6b) inclined at an angle @ against the x-axis, € is a positive quantity 
which increases with the thickness of the airfoil. € is zero for a flat plate and 
is infinite for a circle. Furthermore 


a=i+e. (36) 


The origin of the z-plane lies at the center of the airfoil (Fig. 6b). If we repre- 
sent the circle in the ¢-plane by 


f= ae”, (37) 


the trailing edge of the airfoil corresponds to the point where 6 =0. The lead- 
ing edge of the airfoil corresponds to the point where 6=7. Therefore, the 


chord is 
e 
= 4(1+ p 38 
. ( a+ ;) (38) 


e/a =X 





With 


the values of x and y corresponding to an arbitrary value of @ are given by 
Eq. (35) as 


f cos (@+a)—2Acosa@ 2h? cos “\ 











x = a<cos (0 — a) — A\cosa+ - 
\ a1 —2dcos6+22). 142 i 
. sin(@+a)—Asinéd 2d sina 
y = assin (@— a) + Asina — —_— — ——_—_ 
a?(1 — 2d cos 6 + i?) i+,” 


After Eq. (7), the velocity due to YW» at a point in the z-plane corresponding to 
¢=ae" is 


’ : Ka ( sin (@+ a) —Asina Asin a 
U = v,| 1 + ——}sin (@— a) — — — + \ (40) 
c a2(1—2rcos8+*) a(1 +A) 





This velocity is horizontal and has a tangential component in the counter- 
clockwise direction equal to 


dx 
_ u| : | (41) 
V (dx)? + (dy)? paacl? 
and a normal component, directed along the outward normal, equal to 
dy 
uo ——— =| é (42) 
V (dx)? + (dy)*Jr—aci® 


On the other hand if y; is the disturbance stream function, the velocity com- 
ponent normal to the circle ¢ =ae* is 
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oe s 


Therefore, the corresponding velocity component in the z-plane, normal to 
the surface of the airfoil, is 


, BT 
Then the boundary condition at the airfoil surface requires that this normal 
component of the disturbance velocity be equal to the negative of the normal 


velocity component due to Wo. By means of Eqs. (42) and (44), this boundary 
condition is expressed by 


dg 


rs (44) 








rmae'? 

















1 oy, dt | dy 
C2) Sl--laetele 
r 80 J rma! d2| paaci® V (dx)? + (dy)*Jrmaci? 
But 
. adé 
[V(dx)? + (dy)*];-ce@ = —————- - (46) 
dg 
dz t=ae® 








With the aid of Eqs. (39), (40), and (42), Eq. (45) can be re-written as follows 


1 (- ~*) 
Cert Ose 
ens [1+ "4 sin (@ — a) = alt an ccetind Asin @ \] 


c \ a*(1 — 2rcos8+*) a(1 + dr) 
cos (@ + a) + A cos (@ — a) + 2A cosa 
{cos (@—a | 





(47) 





a?(1 — 2d cos 6 + d*)? 


This equation alone would not determine the function y; completely but for 
the additional so-called Kutta-Joukowsky condition, which fixes the strength 
of the circulation over the airfoil. 

6. Strength of circulation. The Kutta-Joukowsky condition states that 
the velocity at the trailing edge of the airfoil must be finite. The velocity at 
the trailing edge of the airfoil consists of two parts: one part is that due to Po 
and the other that due to ¥:. Only tangential components need be considered 
because the normal components cancel each other as required by the bound- 
ary condition. The part due to y; in the counterclockwise direction is 


_ ah dt 
Or | dz 


(48) 








f=a 


Therefore, by means of Eq. (41), the resultant velocity at the trailing edge is 
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OW dx 1 | dg 
_-—_— — a —_— ; (49) 
or V/ (dx)? + (dy)? | dt dz\¢a0 
dz f=—a 








but {=a is the singular point of the transformation from ¢ to z as can be 
easily verified. In other words, | dt /dz| tends to infinity at =a. Therefore, 
the resultant velocity at the trailing edge can only be finite if the quantity 
within the square bracket of Eq. (49) vanishes. By means of Eqs. (39), (40), 
(46), this condition can be written as 


" (~) — vo] 1+" fsin ee sin (@+a)—A cos t. A sin a \] 
os C 











or a*(1—2d cos 0+?) = a(1+A) 
=0 
sin (@+a)—? sin (@—a)—2X sin a@ 
| sin (@—a)+ — ] =0. (50) 
a?(1—2X cos 0+?) m0 
The appropriate general solution y; in the {-plane is 
. 1 
¥i = Up | a log r + Bo + >> (an cos nO + B, sin n6) =|. (51) 
n=] ig 
Therefore, 
1/1 dy B 2 
a{- - ) = — + >) (= na, sin nO + nB, cos n6) , 
Uo r 00 r=a a n=1 g*tt 
(52) 
1 /oy a = 
—_— (=) ee + >> (na, cos nb + nB,, sin nb) — - 
Uo \ Or J rma a n=l ar 


By expanding the right hand side of Eq. (47) into a trigonometric series, 
all the coefficients a, a2, - + + and Bo, Bi, - - - in yy can be determined by the 
first of Eqs. (52). Then the second equation of Eqs. (52) together with Eq. (50) 
will determine the value of a. Actually, calculations are easier if each 
term on the right of Eq. (47) is taken separately. For example, the term 
— U, cos (@—a) can be expanded into — Uy) cos a cos 8— Up sin a sin 6. Then 
according to Eq. (52), the contribution to a; is a? sin a; the contribution to p; 
is —a*cosa. Finally the contribution to —(0y/0r),-. is Uo(sin a cos 6 
—cos a sin #) = — Ug sin (@—a). The other more complicated terms can be 
treated in a similar manner with the aid of the following expansions and their 
derivatives with respect to @ and X: 

sin 0 2 
> A"! sin n0 
1 — 2A cos 8 + 2 ah 


cos@—X - 
— — = }°r*' cos nb 
1 — 2X cos 6 + i? Feat 4 


1 1 - 
1+ 2 n» cos n?). 
alt +22 


n=1 


II 











1 — 2X cos 6 + i? 
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The value of — (0y,/0r),.. thus found is 


1 Oy ao Kx 
ECD) pt fais 
Te \ Or Sean a c(1+A) 


_ _ sin (@+a) —d? sin (@—a@) — 2, sin “| 
-| sin (@—a) — 
7 @) a?(1—2d cos 6+?) 





Ka 
— E cos 2(0—a) +4 sin (@—a)— 





sin (@+a)—A cos a 
c a?(1—2d cos at 
sin (@+a)—? sin (@—a)—2A sin a@ 
~— @®(1—2d cos 0-22)? 
cos @- ~ a) ) [cos (0+a)—X cos a| 1 
a*(1—2d cos 6+?) 2a*(1—?*)(1—2A cos TH (3) 











The strength of circulation ap necessary to satisfy the Kutta-Joukowsky 
condition can then be calculated by means of Eq. (50). The result ‘is 


n [2 1 = 2X sin? a@ — . 1 ~ +5 1 \] (54) 
= aU | 2s — co -—+—— 
a =a in a 1+) s 2a mad +») 


The strength of the source, Bo, is zero, as would be expected from the fact 
that the airfoil is represented by a closed contour. 

7. Strengths of doublets and quadruplets. By collecting the contributions 
of the different terms of [(1/r) dy1/00],-. in Eq. (47) to a1, ae, 8: and Be, the 
following values of the strength of doublets and quadruplets in the ¢-plane 
are obtained 


1 
au=0'V| (1+) sin a 
a? 
Phen , 4 +(1+ =) + | 
apes = (COS a= enemas endl beneereeenren 
¢ \2e® Ss XA —-2*) a?) a(1+2) 


a 1 Ka . » 1 » 
a=atvsl CH) cos a-}—— sin “at (“Jaa 


rN {- 1 Xe 1—3d? 
ao=a*Uo | sin oye cos 2a ais +— 














a? C 2a? 2a? 4a*(1—d?) 
a & Ke ee ( a \) om 
=a*Uy| — cos a+— sin ee arn neem . 
ailtietied: * oamuuadl™ “V4 2a? 4a(1—2?) 


However, to calculate the lift and moment over the airfoil, it is not the 
strength of doublets and quadruplets in the {-plane that is needed but the 
strength of those in the z-plane, where the airfoil is located. From the known 
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values of a, - - - 8 given by Eq. (35) and the transformation specified by 
Eq. (55), the desired quantities can be easily calculated in the following man- 





ner. 
By introducing the conjugate function ¢; of the stream function one can 


write 


(56) 





bi thy = legend des “= 1a} a hid 12 
e? 


With £=re", the real part of Eq. (56) is Eq. (51) since Bo is found to be zero. 
Now Eq. (35) gives 


' E + : + ] (57) 
t= seia sin eh aman Te 0's 57 
(1 of d)zet# gretia 


for sufficiently large values of z. Then ¢:+7¥1 can be expanded into a series 
in z by substituting Eq. (57) into Eq. (56). The result is 





+i: = ia | +| ati( +=) “tenia 
dit Wi= 1a log z it? ay iA a *é 


+[ b+ Bi 4 ( ay da )| tg-tia (58) 
—B2+——+i[ a2 ———— ap -———_) | x*e*'ee -- - é 
“wa. . ik”6 Uhhh 


If the strength of the circulation, the doublets, and the quadruplets in the 
z-plane are denoted by do, a, 51, a2 and be, respectively, the following relations 
are obtained by comparing Eq. (58) with Eq. (56): 





do = a, bo = 0 
Aa 
1+. 





a= Ai sina + (ai + 
(59) 





1+, 


Bi ; hay a9 
a= (6: — +) sin 2a + as = ay — 08 2a 
+2 1+, 2(1 +r)? 


Aa 
bs = Bi.c0s a — (ai + ) sin a 





Bi Aa Na . 
bs = | Bs — —— } CoS 2a — | ag — —— — a — —— ) sin 2a 
1+A) 1+’ 2(1 + dA)? 


Eqs. (54), (55), and (59) give all the necessary data to calculate the forces 
and moment acting on the airfoil. 

8. Lift and moment coefficients. By means of Eqs. (32), (33), and (34), 
the drag, lift, and moment about the origin can be calculated with the values 
of do, a1, G2, Do, b;, bg given by Eq. (59). It is seen that the drag is also zero in 
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the case of shear flow. The result of the computation can be conveniently rep- 
resented by defining the functions lo, 1;, le, 13, 14 and mo, m,, m2, ms, m4 in con- 
nection with the lift coefficient C, and moment coefficient Cy, in the following 


manner: 


y 
CL= a 2x |lo sin a+K {l:+/2 cos 2a} 4+K? {ls sin a+, sin 3a} | 
2PU oC 
—_ Mo 7] mMo , 
Ca, 2° = sin 2a + K {m, cos a+m; cos 3a} (60) 
3pU,c? 2 2 


+K?} mz sin a+m, sin ta} |. 


The negative sign in the moment coefficient is introduced in accordance with 
the usual convention of taking the stalling moment as positive. The functions 
lo, li, le, 1s, 1g and mo, my, m2, m3, mM, are given by the following equations: 


























1 1 r2 
hb=—) h=— 
h a 1+, 
, - MAt 2/a) (2 + 1/a) 
i= at, 
8h?(1 + Xd) 8h?(1 r 
2°7( n?(1 + dr) 61) 
= — ——— Fes - a+r] 
64h?a(1 + A) L a2(1 + A) 
d 
a ane 
64a2h? (1 +)? } 
1 1 
Mo = - | - | 
ah? a 1+ A) 
1 1 r r 
m, = | —— ~(14a4. —)] 
Sah®La (1 + X) a(1 +d) 
1 1 N 2 2+1 
m= — | (1— 3) + — (5 — 2 *")] 
Sah* La 2(1 + A) \a? 1+ A 
> (62) 
1 | 1 ( 1 ? ) ? | 
m3;>=> — —{—-+ — ee ee 
64ah*L1 + Xd \a? 2a7(1 + A)? 2a(1 + A)? 
. | ( ) A*(1 + 2?) 
mM, = - 1+ vA -— ——_ | 1+ 
64ah! : a2(1 +) a(1 + 2) 2a(1 + r)(1 — 22) 
1 ( 1 mi 1 )+ 9? ] 
4a*\a 1+2% 4a%(1 + r)] ) 
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Table I gives the numerical values of these functions together with the 
thickness ratio 6 of the airfoil. It is seen that with increasing thickness the 
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Fic. 7. Ratio of the distance d from leading edge to aerodynamic center (positive when it 
is behind the leading edge) to the chord ¢ as a function of thickness ratio 6 and non-dimensional 


velocity gradient K. 
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effect of the velocity gradient K becomes larger and larger. For example, if 
5=11.79% the effect of K on C, at small angles of attack is approximately 
given by the term 27K(1,+/,.) =27K X0.02695. In other words, it is equiva- 
lent to a shift in the angle of zero lift by 1.54 K degrees. For an airfoil of 
21.50% thickness, this value is 3.21 K degrees. 





i ad 


Fic. 8. Moment coefficient Cy, , about the aerodynamic center (positive for stalling 
moment) as a function of thickness ratio and non-dimensional velocity gradient K. 


9. Aerodynamic center. To demonstrate the effect of the velocity gradient 
on the moment more clearly, the aerodynamic center for small angles of at- 
tack will be calculated presently. For small values of a the expressions for Cz, 
and Cy, of Eq. (60) can be simplified to 


Cr = nl {lo + K%(ls + 3) }a+ K(, +h)] 


rs : (63) 
Cu, = = [{ mo + K2(2ms + 4m,)}a + K(m; + m:)]}. 


The moment coefficient corresponds to the stalling moment about the 
center of the airfoil. If the moment is referred to a point on the chord at a 








148 HSUE-SHEN TSIEN 


distance d back of the leading edge, the corresponding moment coefficient Cw 


will be 
as ’ ; 1 d 
Cu = Cr, — Cx(— - =) 
2 c 


kad }mo + K2(2mz3 + 4m4)} a + K(m, + mz) (64) 
j 


1 d 

- (> - <) 2a[ {lo + K2(ls + 3ly) fa + K(i + h2)]. 
2 c 

At the aerodynamic center, the corresponding moment coefficient should be 

independent of the angle of attack a. From Eq. (64), this condition gives the 

distance d of aerodynamic center back of the leading edge as 


d 


1 
alias 


1 mo + K?(2m3 + 4m,) 
4 





(65) 


ae K*(l3 + 314) ‘ 


Cc 
Then the moment about the aerodynamic center is given by the following 
coefficient: 
aK mo + K2(2m3 + 4ma4) 
CMac. = - | (ms + ms) — ——_————__— (1 + 4) I. (66) 
2 lo + K*(ls + 314) 
The numerical values of d/c and Cwy,,. calculated from Eqs. (65) and (66) 
are plotted in Figs. 7 and 8 against the thickness ratio 6 for different values 
of K. 
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ON THE DISCONTINUOUS FLOW AROUND AN 
AIRFOIL WITH FLAP* 


BY 
P. W. KETCHUM 


University of Illinois** 


Summary. A method is given for calculating the effect of the wake on the 
pressure distribution produced by the steady, two dimensional flow around 
an airfoil with infinite span and broken line profile. To determine the size 
of the wake there is proposed the hypothesis that the flow is such as to create 
a minimum disturbance at distant points. 

1. Introduction. The general Navier-Stokes equations for the steady, vis- 
cous flow of an incompressible fluid around an immersed body have not been 
solved even for such simple cases as the sphere or cylinder. However, it is 
well known that certain aspects of the flow are adequately represented by 
neglecting viscosity entirely. Thus the perfect fluid theory of Kutta and 
Joukowski for calculating the reaction ef the air on an airfoil with infinite 
span, engaged in a uniform rectilinear translation through a bulk of air at 
rest, gives satisfactory results so far as the lift is concerned. This theory fails 
to account for all the details of the flow, however. In particular, it neglects 
the wake of deadwater which exists at the trailing edge. The effect of this 
wake is to make the lift slightly smaller than the calculated value and to 
reduce the negative pressures along the rear portion of the suction side. Since 
the effect is greatest near the trailing edge, the discrepancies between calcula- 
tions and observation will be most pronounced for the total hinge moment 
of the pressures on a trailing edge flap. 

The boundary layer theory has been used to calculate the size of the wake, 
i.e. the point where the wake detaches itself from the body. Basing his cal- 
culations on the’ observed pressure distribution, Hiemenz in 1911 obtained 
from the boundary layer theory an angle of 82° from the forward stagnation 
point for the point of detachment in flow around an infinite circular cylinder, 
a value in good agreement with his experiments. However, his experiments 
were made at low Reynolds numbers; an airfoil would be stalled for such a 
flow. For the high Reynolds numbers corresponding to the normal, unstalled 
flight of an airfoil, the calculated point of detachment for flow around the 
cylinder is at least 20° smaller than the observed value, even though the cal- 
culations are based on experimental pressure distributions. 


* Received Dec. 30, 1942. 

** This investigation was undertaken at the suggestion of Professor W. Prager while the 
author was a Fellow under the Program of Advanced Instruction and Research in Mechanics 
at Brown University, 
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A definite advance was made by C. Schmieden in 1932! when he showed 
that for high Reynolds numbers the perfect fluid theory is adequate to ex- 
plain the wake behind a cylinder if one assumes a discontinuous flow. ‘Accord- 
ing to this view, the wake is obtained as a special case of the classical theory 
of jets. The theory of jets, due to Helmholtz and Levi-Civita, assumes a 
wake of stationary fluid at uniform pressure, bounded by two streamlines 
along which the speed is constant and different from zero. In spite of the fact 
that such a discontinuity in velocities would be impossible in a real fluid, 
Schmieden’s calculated values of the pressures around the cylinder agree re- 
markably well with experiment. The agreement is perfect up to the pressure 
minimum; beyond that there are some deviations, which can be attributed to 
turbulence in the boundary layer and circulation of fluid in the wake. The 
point of detachment is approximately correct, but the observed point is some- 
what masked by the turbulence; and the circulation of the fluid in the wake 
causes the constant pressure in the wake to be negative instead of zero as 
predicted by the theory. 

Schmieden found, however, that the equations of motion for a perfect fluid 
plus the usual boundary condition that the velocity approach the stream 
velocity as the distance from the cylinder increases without bound are insuffh- 
cient to determine the discontinuous flow uniquely. Various sizes and types 
of wakes are possible. There exist symmetric discontinuous flows with wakes 
detaching from the cylinder at every angle from about 55° to 180°. For angles 
between 55° and about 124° the wake extends to infinity and is bounded by 
streamlines which are first concave toward the wake and then, after passing 
through a point of inflection, are convex; the drag on the cylinder is positive; 
and the constant pressure in the wake is zero. For angles between 124° and 
180° the wake is finite in extent, the drag is zero, and the pressure in the 
wake is positive.” For the critical angle of about 55° the wake boundaries 
are always convex (Helmholtz case). For the critical angle of about 124° the 
wake is infinite but has a width which tends asymptotically to zero, the 
boundaries of the wake are concave, the drag vanishes, and the pressure in 
the wake is zero; we shall refer to this as the Schmieden case. 

Since a correctly set physical problem should have a unique answer even 
for a perfect fluid, an additional boundary condition is evidently required 
to determine the flow. Schmieden has chosen the auxiliary condition that 
the wake boundaries be concave and infinite in extent. A possible reason for 
selecting the condition of concavity is the intuitive feeling that in a viscous 
fluid any portion of a streamline with a point of inflection would be swept 


1 Schmieden, C., Uber die Eindeutigkeit der Lisungen in der Theorie der unstetigen Strim- 
ungen, Ingenieur-Archiv 3, pages 356-370 (1932) and 5, 373-375 (1934). 

* Schmieden considers only infinite wakes, but it is obvious that finite wakes can exist: 
see Kolscher, M., Unstetige Strémungen mit endlichem Totwasser, Luftfahrtforschung 17, 154- 
160 (1940). The present paper also restricts attention to infinite wakes. 
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away to infinity. However, for a perfect fluid there would be no such tend- 
ency; moreover, the condition of concavity is an “internal” condition, more 
in the nature of a law of motion than of a boundary condition. The required 
boundary condition must be “external” in character—a condition which can 
be imposed arbitrarily as a constraint on the system. Such an external condi- 
tion could only be imposed by specifying something about the nature of the 
impinging stream. In particular we could require that this stream be as nearly 
uniform as possible, that is, the difference between the velocity and the un- 
disturbed stream velocity shall approach zero as rapidly as possible as one 
recedes from the cylinder in the direction of the stream source. In the present 
paper we show that this assumption of “minimum disturbance at infinity” 
does in fact lead to the Schmieden case, at least for infinite wakes. 

On the experimental side the justification for choosing the Schmieden case 
for an approximation to the flow of an actual fluid at high Reynolds numbers 
lies in the agreement previously mentioned between calculated and observed 
pressure distributions.* Of the various discontinuous flows the Schmieden 
case gives the best agreement for the pressure distribution as a whole, but 
more particularly for the forward part of the cylinder where turbulence is 
absent. As previously mentioned the observed pressure in the wake is nega- 
tive instead of zero, but this is still the best agreement possible since the 
theory predicts a positive pressure for finite wakes. Quantative technical ap- 
plications involving pressures behind the cylinder are subject to the limita- 
tions imposed by this discrepancy. Further refinements of the theory to 
eliminate this discrepancy by taking into account the turbulence in the wake 
have been suggested by Schmieden (Reference 1, page 364). 

In 1940 Schmieden‘ applied his theory of discontinuous flow to the case of 
an airfoil with infinite span whose profile is an inclined straight line segment. It 
was found that the pressure distributions so obtained agree almost exactly 
with the pressure calculated from the Kutta-Joukowski theory except near the 
trailing edge, where the wake would be expected to exert an appreciable in- 
fluence. It should be observed that the Kutta-Joukowski theory also uses 
an additional assumption beyond those involved in the perfect fluid equa- 
tions, namely, the Joukowski hypothesis that the circulation is such that 
the velocity at the trailing edge is finite. Schmieden makes no direct assump- 
tion about the circulation, but replaces the Joukowski hypothesis by the 
above assumption concerning the nature of the wake. From Schmieden's 
point of view the Joukowski hypothesis is a special assumption concerning the 
nature of the wake in the limiting case when the wake is made to disappear. 


3 Schmieden compares his results with the experimental values of Eisner, F., Widerstands- 
messungen an umstrimten Zylindern, Mitteilungen der Preussischen Versuchsanstalt fiir Wasser- 
bau and Schiffbau, Heft 4, Berlin 1929. 

‘Schmieden, C. Uber Tragfligelstrémungen mit Wirbelablisung, Luftfahrtforschung 17, 
37-44 (1940). 
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An objection has been raised to Schmieden’s theory of the airfoil with 
straight line profile in that the pressure is negatively infinite at the leading 
edge; it is one of the requirements of the classical theory of jets that separa- 
tion of the fluid from the body occurs at such points. However, in an actual 
airfoil no separation occurs at the leading edge (unless the airfoil is stalled) 
because the leading edge is rounded off. The infinite velocity is not a neces- 
sary feature of the Schmieden theory, but only results from the sharp leading 
edge which was assumed to render the mathematics feasible. Singularities of 
exactly the same order are obtained from the Kutta-Joukowski theory when 
applied to sharp leading edges. 

In the present paper we apply the Schmieden theory to the case of two 
dimensional flow around an airfoil with flap, the profile being a broken line. 
Some of the results are compared with the Kutta-Joukowski theory in Figs. 
2 to 5. 

2. The general results. Consider a plate whose cross section is a broken 
straight line (Fig. 1). The angle of attack of the main wing may be denoted 





Fic. 1. 


by a and the angle of attack of the flap by a+. Thus £ is the angle of de- 
pression of the flap. We assume that a is positive and that 8, if negative, is 
greater than —a. The velocity of the undisturbed stream (velocity of flight) 
is denoted by vq. 

There will be a critical streamline which divides at the stagnation point S, 
follows both sides of the airfoil to the points of detachment, D and T, and 
then to infinity as boundaries of the wake. At the points of detachment the 
boundaries of the wake leave the plate tangentially. 

Evidently the ratio of the lengths (chords) KT to LK must not be too 
small, or else the point of detachment D will not fall on the flap. We restrict 
attention to cases where the flap is long enough for D to lie on the flap—an 
assumption that appears to be justified in all cases of technical interest. 

The method of Levi-Civita* does not give the velocity distribution along 
the surface of the plate directly in terms of the position which the point 
under consideration has on the plate; but instead, both position and velocity 


* For an exposition of the Levi-Civita theory in English see: Milne-Thomson, Theoretical 


Hydrodynamics, London, 1938, pp. 305-316. 
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are expressed parametrically in terms of an auxiliary variable ¢, where ¢ 
is an angle ranging from 0 to m7. The explicit expressions are as follows, where 
v is the speed of the fluid and s is the distance measured along the broken 
line from T. (For the proofs of these formulas see section 5): 


sin $(@,+¢@) sin in 3 | os—o | Ee 3(¢x+¢@) sin 4 } | or— oT - 





sin }| ¢.—| sin }(¢s+¢) Lsin | ¢x—¢| sin }(¢2+¢) 
, * sin 2(d1—¢) sin? } 2(¢s+¢) sin 3| ox— ¢| sin 2(¢e+9) Ale 
1c? ‘ sin gd. (2) 


9) 9 
t Veo 





S=4 ———_————_—___-—— —$—$$$_____—_—__ 


sin }(6.+¢) sin }(¢x+¢) sin }| ¢z—¢| 


Here the constants dz, ds, 6x, @x and c? depend on the dimensions of the 
plate and on @ and 8. The number c? is a proportionality factor |determining 
the scale of the figure. Values of @ between 0 and ¢, correspond to points 
on the under side of the plate while values between ¢, and 7 give points on 
the upper side. The determination of the numbers ¢z, $s, $x, and ¢x is the 
most difficult part of the whole procedure and will be discussed later. For 
the present we suppose that they are known. It will be sufficient now to ob- 
serve that these four constants are angles between 0 and 7, whose magnitudes 
are in the order written, which if used as the upper limit in (2) will yield the 
positions of the points E, S, L, and K respectively. Also, the angles ¢s and ¢, 
are near 1/2 if a and B are small. 

If the integration in (2) is performed there is obtained with (1) a com- 
plete parametric representation of the speed at the airfoil surface. From 
Bernoulli's theorem it is then a simple matter to calculate the pressure in 
terms of the parameter @. If the pressure at infinity is p,, the pressure at 


any point will be 


p=—(ta—0) + po (3) 


2 


where p is the density of the fluid. 

The total hydrodynamic forces acting on all or part of the airfoil may be 
found by integration of the pressure. In particular, the total moment of the 
forces acting on the flap taken with respect to the hinge will be 


M 7 1 ) “ if f- l ) = do (4) 
l = p= § —_— _ p= $ —_— ’ 
. 0 PS . : Phe oK P's do 


where lp is the length (chord) of the flap ET and moments in the counter- 
clockwise direction are considered positive. 
The resultant force on the airfoil is vertical; the drag vanishes and the 


total lift is (section 9): 


9 


22 esas ; 
Y = — 3rpv,..c sin 2¢1 — sin 26s + — (sin 26x — sin 26e) | 
T ‘ 
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In addition to the distributed pressures there is a concentrated force 

(suction) at the leading edge which acts parallel to the plate LK (see section 

10). Its horizontal component must equal the horizontal component of the 

resultant pressure on the airfoil, since there is no concentrated force at the 

hinge K. The magnitude of the force at L is 

sin 3(¢1 + $x) sin 3(@z — sey 


F2x4 ove * sin’ ¢ sin 3(¢ 6s)| 
1 = 47p0.C€ L aes ~ O4)1-;.. : 
sin 4(¢x — $x) sin }(o1 + ¢e) 





The integrations in (2) and (4) must in general be carried out numerically 
or graphically. The integral in (2) has been carried out in closed form only 
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for the case of a straight plate, i.e., 8=0, (Schmieden, Reference 4). The 
problem is much simplified in that case by the fact that of the four constants, 
oe and ¢x do not occur and ¢,=}(47+a), ¢s=}(4—a). (See section 4 for 
the proof). 

3. Some numerical results. If the angle of depression of the flap is very 
small in comparison with the angle of attack, we consider as a first approxima- 
tion that 8 =0. Schmieden (Reference 4) has given a diagram for the pressure 
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Fic. 3. 
Lf =Total lift on flap without wake, lp=length of flap= KT; ae 
Lr” =Total lift on flap with wake, 1 =total length of airfoll=LA+AT. 


distribution in this case for an angle of attack of 12°, showing a comparison 
of values obtained both with (Schmieden theory) and without (Kutta- 
Joukowski theory) a wake. His diagram is reproduced in Fig. 2, together with 
a similar diagram for an angle of attack of 3°. 

Fig. 3 shows the effect of the wake on the total lift Ly of the flap, for 
flaps of various sizes. Fig. 4 is a similar diagram for the hinge moment My 
of the flap. 

Fig. 5 shows the effect of the wake on the pressure distribution for the 
case a=10°, 8=5° and \=KT7/LK =1.00.5 


5 For the Kutta-Joukowski solution see: Schmieden, C., Die Sirémung um einen ebenen 
Tragfliigel mit Querruder, Zeitschrift f. angew. Math. u. Mech. 16, 194-198 (1936). 

See also Ellenberger, G., Bestimmung der Luftkrifte auf einen ebenen Tragfliugel mit Quer- 
ruder, Zeitschrift f. angew. Math. u. Mech. 16, 199-226 (1936). 
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4. The conditions for the determination of the constants. The Schmieden 
hypothesis. As ¢ varies from 0 to dg the corresponding point in Fig. 1 will 
move along the under side of the flap from T to E; as ¢ goes in succession 
from ¢z through ¢s, ¢z, ¢x to 7, the corresponding point will move along 
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My =Hinge moment on flap without wake, Jy=length of flap=KT, a 
My" =Hinge moment on flap with wake, 1 =total length of airfoll=LA+AT. 


the under side of the plate through S to Z and then back along the upper side 
through K to D. Now the distance traversed on the under side between E 
and LZ must equal the distance on the upper side from L to K. Hence, if the 
integrand in equation (2) is f(¢, dz, s, dz, Ox), then 


>] 


PL PL 
f S(o, de, os, Ox, Ox)dd = f Sd. (5) 
OK 


E 
If \=7K/KL is the ratio of the length of the flap to the length of the station- 
ary part of the airfoil, we must have 


i fdp = de. (6) 


OE 


Equations (5) and (6) are the two conditions on the constants imposed 
by the geometry of the airfoil. A third condition is obtained from the require- 
ment that the velocity must approach the stream velocity as the distance 
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from the airfoil tends to infinity. The following equation is the mathematical 
formulation of this condition (for the proof see section 6): 


B 
o. — ¢s + — (ox — oe) =at+B. (7) 
T 


For the fourth condition we need the Schmieden hypothesis mentioned 
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in the introduction, or some alternative. The hypothesis which we propose 
is a natural and direct extension of the physical assumption on the basis of 
which (7) was obtained. Whatever physical principle causes the velocity to 
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approach the stream velocity as the distance from the airfoil becomes large 
might also be expected to cause it to approach the stream velocity as rapidly 
as possible. This is our basic hypothesis; it says that the disturbance at distant 
points shall be a minimum. 

The general mathematical formulation of our hypothesis will be found in 
section 7, and in section 8 it is shown that for our case of a broken straight 
line the hypothesis leads to the condition 


sin @, — sings + La (sin dx — sin dz) = 0. (8) 
us 


For 6=0 this reduces to the condition which Schmieden (Reference 4) has 
used for determining the size of the wake, namely ¢,=7—@s. Putting 6B =0 
1 


in (7) we get 6. —¢?s =a. Hence ¢,=}3(x+a) and ¢s =3(4r—a). 

The four equations (5), (6), (7), and (8) will determine the four constants 
involved: These equations can only be solved by successive approximation. 
Slide rules can be constructed for solving (7) and (8). Equations (5) and (6) 
can be solved by trial, the integrations being performed either graphically or 


numerically. 
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5. Proof of formulas (1) and (2) for the velocity distribution. We consider 
Fig. 1 as a diagram in a complex z-plane, z=x-+2zy. We choose the complex 
potential 
w=~it we 
in such-a way that the streamline ¥2=0 will be the critical streamline which 
follows the surface of the airfoil and wake, with yi1=0 at S. 

The w plane minus the positive half of the real axis is mapped onto the 
interior of the upper half of the unit circle in a 7 plane (Fig. 6) by the trans- 
formation (for details see, for instance, Reference 4): 


(+4) ° 


Le 


w= a* (cos os — - 
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The radii from —1 to 0 and from 0 to 1 in this r-plane correspond to the free 

streamlines. The points on the unit circle with arguments @¢z, ¢s, $x, and ox 

correspond respectively to E, S, ZL, K. The number a is a constant. 
Consider the complex velocity 


dw 


—— = os, 10 
dz (10) 


Along the real axis in the r-plane between —1 and +1 the speed v must be 
constant and equal to v,,, the undisturbed stream velocity. Along the succes- 
sive segments of the unit circle between the points 7’, E’, S’, L’, K’, and D’, 
the angle 6 must have the respective constant values —(a+ 8), —a,7—a, —a@ 
and —(a+ 8). It will now be shown that the following expression for the com- 
plex velocity will satisfy these conditions: 


dw OM, eh: ‘eu)eten(r — [F — e**x)e'*x(r — |" (11) 


(r — et*x)(r — e~ **z) et*e 





dz (r — e*1)(7r — ec **s)e**s 


We use that determination of the bracketed quantity to the 8/7 power whose 
argument vanishes for tr = 1. The absolute value of this expression is obviously 
constant and equal to »,, if 7 is real. For r =e on the unit circle we make use 
of. the identity 


ei — eit’ — eli(ot+e’) (ehil(e—o’) on e~ ti(e-4")) 


Ziedi(et+o” sin 3(¢ — ¢’) 
to rewrite our expression in the form 


dw i sin 3(¢ + 2) sin 3(¢ — $s) 
= v,,€°'* _ _—$ $$$ $$ 


dz sin 3(@ — 1) sin 3(@ + $s) 


{= 3(@ + x) sin 2(¢ - =< i iii 


sin }(@ — $x) sin 3(¢ + ¢z) 


where the determination is chosen to agree with that in (11). If ¢ is between 
0 and ¢ each factor involving the sine of a difference is negative, all the nega- 
tive signs will cancel and the argument of the whole expression will be a+ 8; 
hence 6= —(a+8) as desired. If ¢ is between ¢g and ¢s, the quantity inside 
the brackets will be negative and have an argument —7; hence the entire 
expression will have an argument a. In this way all the prescribed values of @ 
may be verified. Since our expression for dw/dz is analytic and not zero in- 
side the upper half of the unit circle in the 7-plane and satisfies the required 
conditions for v and 6 on the boundary of that semicircle, it follows that the 
flow around our airfoil is thereby correctly determined. If we take absolute 
values of both sides of (12) we get the required formula (1). Also, from (9), 


for r=e*, 
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dw 
— = — a*(cos ds — 3(7 + 1/7))(1 — 1/7?) 


dr 


Il 


— 2a?i(cos ds — cos d)e~** sin d 


4a*ie~** sin (ds + ¢) sin 3(¢s — ¢) sin ¢. (13) 


Il 


Since dz = (dz/dw) (dw/dr) dr, we have 
; 1 | dw |, 
| dz| = —|—|| do]. (14) 
v | dr | 
Now ds will have the same sign as d@ if ¢ is less than ¢, and opposite signs 
if 6>dz. We can therefore get ds from | dz| by removing the absolute value 
signs around d@ and around the quantity ¢,—¢ in the expression (1) for v. 
Hence (2) follows immediately on substituting (1) and (13) into (14), in- 
tegrating, and setting c?=a?/v,.. 

6. Proof of equation (7). We need only observe that the complex velocity 
must approach v,, as T—>0. Putting 7 =0 in (11) we get 


Vx = U,€*s eth € -i¢zeits[e-i*xeiee |B ¥. 


which reduces at once to (7). 
7. Formulation of the hypothesis of minimum disturbance at distant 


points. Consider the function 


dz Veo 
w(r) = In (>. ) = In— + #. (15) 
au 


e 
v 


This function is regular within the upper half of the unit circle and continu- 
ous on the diameter formed by the segment of the real axis between —1 and 
+1. It is also pure imaginary on this diameter. Hence it can be continued 
into the lower half of the unit circle by reflection of iw(r); and w(r) will then 
be regular within the entire unit circle. 

Expanding w(r) in a power series we get 


dw F ° 
— w(t) «=. ¢ -w (0)—w’ (0) r—w" (0) r*/2—- +> 
——— Vao€ —_ Vin€ w(0)—w’ (0)7 ))r*/ 


dz 


7? 
- 24 — w'(0)r + [(w’(0))? — w"(0)] iene \. (16) 


Here we have put w(0) =0, which is the condition that the velocity equals v.,, 


at infinity. 
In all cases discussed in detail by Schmieden and in the present paper, 
the flow is uniquely determined by the condition 


w'(0) = 0. (17) 
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Hence, there is one and just one flow for which 
-—G~r of ¢~~G 


for n22; for every other flow n=1. Since r-0 as z— (for points outside 
the wake), dw/dz will approach v,, at a maximum rate as z—© if the flow 
satisfies (17). Thus (17) is, for the cases under consideration, the desired 
mathematical formulation of our hypothesis that the disturbance of the main 
stream at infinity shall be a minimum. 

Equation (17) would be insufficient to determine a unique flow if one of 
the two points of detachment of the wake had not been determined a priort. 
Thus in the case of an airfoil with sharp trailing edge, one point of detach- 
ment is placed at the trailing edge; for the case of a circular cylinder, the two 
points of detachment are assumed to be symmetrically placed. Situations can 
be imagined, however, such as an airfoil with blunt trailing edge, where 
neither point of detachment could be determined in advance. In such a case 
one can have a still smaller disturbance of the main stream at infinity (i.e. 
n=3 in the above asymptotic formula) by imposing in addition to (17) the 
condition 

w”(0) = 0. (17’) 


For symmetric flows this last condition is automatically satisfied whenever 
(17) holds. 

The theory predicts the following sequence of events as one passes con- 
tinuously from the case of a circular cylinder to an inclined plate by way of 
elliptic cylinders with increasing eccentricity but with major axis having a 
small constant inclination with the direction of the stream: for small eccen- 
tricities (and very high Reynolds numbers) both w’(0) and w”(0) vanish, both 
the drag and lift are zero (see section 9), and the disturbance of the main 
stream at distant points is the least possible. As the ellipse becomes flatter 
and the trailing edge less blunt, the local viscosity and pressure conditions 
at the trailing edge gain control over the lower point of detachment of the 
wake, the condition (17’) is lost, and a corresponding lift is developed. The 
flow in this state corresponds to that in the normal flight of an airfoil. For 
still flatter ellipses or higher angles of attack or lower Reynolds numbers, the 
local conditions around the leading edge gain control of the upper point of 
detachment, the condition (17) is also lost, and there is a positive drag propor- 
tional to w’(0) and a positive, although much smaller, lift. In this state the 
airfoil is clearly stalled. Thus the development of a lift and of a drag mark 
abrupt changes in the state of the flow associated with the loss of the condi- 
tions (17’) and (17) respectively; and there is a critical range of Reynolds 
numbers marking the effect of local conditions which determine the transi- 
tions from one type of flow to another, 
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8. Proof of equation (8). Consider the expression 
wi(r) = In (r — e*#1) — In (7 — e-*#1), 


where ¢; is an arbitrary constant angle. Then 





= 1 1 
wi (r) = - —_— 
T= e'?i y eel 
Wf 1 1 
w"(r) = — ——eemeienans Te 
(r ane eter)? (r = e~i¢1)2 
wi (0) = — e~*#1 + eft: = 27 sin gu, 
wi"(0) = — e261 + e?t61 = 27 sin 291. 


Now from (11) it follows that w(r) is the sum of a constant and constant 
multiples of expressions of the form w:(r). Hence 


w'(0) = 23 sin oz — sings + £. (sin ¢x — sin és) |, 
¥ (18) 


: B 
w"(0) = 2i| sin 21 — sin 2¢s + — (sin 26x — sin 2¢z) |. 
Tv 


Equation (8) is an immediate consequence of the first of these equations and 


our hypothesis (17). 
9. The total drag and lift. According to the Levi-Civita formula, the total 


drag on the airfoil is proportional to w’(0). From our hypothesis (17) it fol- 
lows that the drag is zero. In fact, in the case under consideration our hy- 
pothesis is equivalent to the condition that the drag vanishes. However, in 
other situations such as the case of a finite wake discussed by Kolscher 


(Reference 2), the two conditions differ. 
The total lift may be calculated from the Levi-Civita formula 


22, 2 2 
VY = inpv,..c iw” (0), c =a/t.z, 
and (18) to be 
‘ oof . ioe ; 
VY = — 4200..¢ sin 261 — sin 26s + — (sin 26x — sin 261) |. 
T 
For 8 =0 this reduces to the Schmieden expression 


4... 
Y = mpv..c sina 
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10. The concentrated forces at ZL and K. Consider the flow around a 
corner K with arbitrary angle yr, 0Su<1. Draw the small circle C of radius r 
about K. The finite force at K, if any, will be obtained by applying Euler’s 
momentum theorem to the fluid inside C and taking the limit as r approaches 








Fic. 7. 


zero. In Fig. 7 we show the fluid as a free body, where — F is the reaction of 


the concentrated force F on K. 


We have 
e-f ipds= pf wads, 
C+AKB ¢ 


where u is the inward unit normal. 
It is assumed that the stagnation point S is not at the corner. Hence 


for small r, setting \=yu/(1+y), OSA SH, 


dw ‘ 
— = k/(e — KP + 0) (19) 
where & is a constant and 0(r'~*) is a quantity such that 0(r!-*)/r!-® is 
bounded as r-0. Hence if we put z—K =re*¥, 
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9 


y=|k| r+ O(r?-), v | &|27-™ + O(r'-), 
ce? = (| k| /k)e®* + o(r'-), ~P=- 3 | k |? pr-? + O(r’), 


where ¢ = 0 if \S1/3, =1—3) if X>1/3, and o(r!—*) is a quantity such that 
o(r!-) /ri-2— 0 as r-0. 


Now 
K ee r x r 
f npds = — 3n| kite { r—dr +n f O(r*) dr. 
{ 0 0 


Both these integrals tend to zero with r if u<1. If u=1, this integral and the 
corresponding integral along KB both diverge, but since, in this case, the 
pressures along AK and KB are directly opposing, it is permissible to use the 
principal value of the integral along AKB. With this understanding we may 


f n pds = — An, | kite f (x! — r“)dr +f n O(r*)dr, 
iKB 0 0 


where ; is the normal to AK. Both these integrals tend to zero. For the 


arc C, 
n ds = — 5 k\*p r—At+lei¥—mqy + n O(r?) ds. (20) 
p ' ' 
v7 0 ¢ 


These integrals also approach zero; the first integral is zero if u=1. 
The Euler momentum theorem thus reduces to 


F = lim ef n-vvds. 
r—0 C 


But cos (n, v)=cos (W—0—7), 1=ve*?, 


= ns = p 9 ] p 9 9 
rf n-vods = — f vei¥ds — - f ve ‘ve? "ds, 
c 2 rf 2 4 


The first of these integrals is similar to (20) and will approach zero for the same 


reason. For the second integral, 


a r(lt+u) | B|2 | k|? 
p m Par p P Pp 
—-= f se ‘¥e2"ds = f e~*¥ —_—_ ¢?)¥rqdy — £ fo )ds. 
L2Jc 2 0 "sash k? 2J7¢ 


write 


The last integral approaches zero. The first integral also approaches zero if 
u<i. If ~=1, however, its value is not zero. Evaluating this integral we get 


finally 
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We have thus proved the general result that for a concentrated force to 
exist at a corner which is not a stagnation point the corner must be a cusp. 
In particular, there will be no concentrated force at the hinge K. The con- 
centrated force at the leading edge may be found by evaluating k, giving 

— sin 4(@, + ox) sin $(¢z — ser 


> ie 2 2 
F = —4pmv,ce sin @, sin 3(¢1 — $s) Ee —— 
_ sin }(¢x — x) sin 3(¢z + ¢z) 





where, as before, c?=a?/v,,. The direction of the force is parallel to the plate 
LK and toward the left. In the case 8 =0 this reduces to 


2 —ta 


9 
F = — prv.c € 


‘ _ 2 
sin a, 
which agrees with the expression given by Schmieden (Reference 4). 

11. The width of the wake. It will now be shown that under the minimum 
hypothesis (17) the width of the wake tends to zero asymptotically as the 
distance from the airfoil becomes infinite. 

We have as in (16) 


dz 1 1 7? 78 
= —¢e"') = — E + w”(0) — sorry |, 
w 2 6 


du Ve Veo 
Also 
dw - 1 1 1 , 
de = — q? 73 CS cos ds + cos ds — _ m 
Hence 
az d u a’ 1 1 1 w” (0) 
dz = dr = —- — . — — cos os + —— — + 0(1) lér 
li nid Ve “ T 7* 


a? sy 7 w"” (0) 
(et, +s | - + eS eee 0 |, 


4 7? T 


where C; and C; are real constants. This expression is of the form 
i 
— + —--+ iD in ¢ + Or) 


TT T 


2=C, + @,+ 


where A, B, and D are real constants and A is positive. For 7 real and positive 
this will give the coordinates of a point (x1, 1) on the lower boundary of the 
wake: 


AB 
x = Ci +—+—+4+0(), 
a a (22) 


y1 = C2 + Din + Or), 
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If r is negative we put r= —7’ where 7’ is positive and obtain in like manner 

the coordinates (x2, yz) of a point on the upper boundary of the wake: 
AB 

2 =C, + eee ee See ae aD + O0(7’), 


T T 
vo =Co2+ Dine’ + O(7’). 
Consider the particular values of rt and r’ defined by the equations 


1 — B+ /B?—4AC,-—) 

T 7 2A 

1 B+ /B? — 44(C, — xD — 8) 
i 24 : 


7 


where ¢ is a large real positive parameter. As ¢ becomes infinite both 7 and r’ 
tend to zero through positive values. Substituting these expressions into (22) 


and (23) we get 


Xe — x, = O(r’) — O(r), 
— B+ /B? — 44; — 0d) 
yo- VN = D In ——— — - + O(r’) om O(r). 
B+ +¥VJB? — 4A(C; — xD — 2) 


As i both these differences approach zero; hence the width of the wake 
tends to zero. 

12. Details of the calculations. In the numerical calculation of pressure 
distributions, the evaluation of the constants is the most time consuming. 
As mentioned previously, it is a simple matter to construct slide rules for 
Eqs. (7) and (8). It was also found to save time to evaluate the integrals 
in (5) and (6) graphically. Templates were constructed for the functions 
log sin ¢, log sin ¢/2, and 8/7 log sin ¢/2. By properly adjusting the phase, 
one can quickly draw from these templates the curves log sin 3(¢—¢z), 
log sin 4(¢+¢ 2), 8/7 log sin 3(¢@ —¢z), etc. By addition and subtraction of the 
ordinates thus obtained one gets the logarithm of the integrand, and hence 
the value of the integrand itself from a graph of the logarithm function. 
The integral can then be found either graphically or by means of a planimeter. 
The pressure distribution shown in Fig. 5 and the associated points in Figs. 
3 and 4 were calculated in this way. 

If the integrals in (5) and (6) are computed numerically, the infinite dis- 
continuities of the integrand at ¢g and ¢x are troublesome. There is not 
actually a discontinuity at ¢x, but the integrand goes to zero so abruptly 
near this point that it is convenient to treat the point as if it were a discon- 
tinuity. The usual method of “subtracting off” the discontinuity is unsatis- 
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factory because the result will be an abrupt zero of the type just mentioned. 
The following procedure was used as an alternative: after isolating the dis- 
continuity (say at ¢,) at the center of a small interval (¢:—6, ¢:+4), and 
approximating the sine of an angle by the angle, the integrand reduces to 
the form g(@)|¢—¢:|’, where v is some constant and g(@) is continuous (and 
smooth). A linear transformation will change this integral into one of the form 
S76") ¢' \*d@’. One can then approximate 2(¢’) by a quadratic function 
of the form ag’*+)¢'+<c by interpolating to the values of g at —1, 0, and +1; 
and the resulting integral is easily computed. The final result is a formula 
for the original integral in terms of the values of g at ¢:—4, gi, and ¢1+6 
which is strictly analogous to, and in fact is a generalization of, Simpson’s 


rule. 
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ON A CLASS OF DIFFERENTIAL EQUATIONS IN 
MECHANICS OF CONTINUA* 


BY 
LIPMAN BERS anp ABE GELBART 
Brown University 
1, Introduction. This paper is concerned with partial differential equa- 
tions of the form 
uc = T1(y)0 
a (1.1) 


Uy = — T2(y¥)V2, 


and with the second order equations 


TWer T (rivy) = 0, 
which are obtained from (1.1) by eliminating either v or u. Here u=u(x, y), 
v=v(x, y), and subscripts x, y denote partial derivatives. Equations of this 


c 


form are frequently found in problems of mechanics of continua (cf. the ex- 


amples in sections 5-7). 
The system (1.1) possesses the same structure as the Cauchy-Riemann 
equations 
Uz = Vy 
(1.3) 
Uy = — Vz 


connecting the real and the imaginary parts of an analytic function of x+7y. 
This similarity suggests an integration theory similar in pattern to that of 
the complex function theory. The fundamentals of such theory are presented 
in this paper. (A more elaborate mathematical treatment, containing all 
proofs, will be published elsewhere). The theory will be illustrated by some 
physical examples. In treating these examples our aim is not to obtain new 
results in mechanics but rather to present known facts from a simpler and 
more unified point of view. 

In what follows we suppose that the coefficients 7; (¢=1, 2) are positive 
analytic functions of the real variable y. Then the equations (1.2) are of 

* Received Feb. 16, 1943. This paper was written while the authors participated in the 
Program of Advanced Instruction and Research in Mechanics at Brown University, Summer 
1942; it was presented to the American Mathematical Society on October 31, 1942 under the 
title, On solutions of the differential equations of gasdynamics. The authors are indebted to Pro- 
fessor W. Prager for his suggestions, criticism and constant encouragement. 
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elliptic type with analytic coefficients. Therefore u and v will be analytic 
functions of the real variables x and y. 

2. =-monogenic functions. If u and v are two conjugate harmonic func- 
tions, (i.e. solutions of (1.3)), then f=u+1v is an analytic function of the com- 
plex variable z=x+iy. From u and v we can obtain two other pairs of con- 
jugate harmonic functions, U, V and w’, v’, by setting 


* dj 
U+ iv -f fdz, u’ + iv’ = f(z) = —> 
20 dz 
i.e. 
U = f udx — vdy, V -f vdx + udy 
z0 20 
and 
u’ = Uz = Vy, vy =9, = — ty. 


This procedure can be extended to equations (1.1). Let u and v be a pair 
of solutions of (1.1). We define the functions 


z z Uu 
U -{ udx — Tovdy, V -f vdx + —dy (2.1) 
z T1 


z 
20 0 


u = Uz = T1Vy, “w=7,=-—: (2.2) 


In (2.1) the integration is extended from an arbitrary fixed point x9+tyo = 2 
to a variable point x +iy =z. By virtue of (1.1) these line integrals do not de- 
pend upon the path, but only upon the points 2, 2p. 


We have 


u 
U,=u, Uy = — 72, Ve-=27, Vy=—»> 
T1 
and 
4 \ ‘ Urry ‘ 
Uz = T1V zy, Uy = Uy, zZ = ’ Vy = Vey 
T2 
Thus U, V and u, v are solutions of (1.1). 
It is convenient to consider 
f=utia 


as a function of the complex variable z=x+iy. If u and v satisfy (1.1) we 
shall call f a 2-monogenic function, = denoting the matrix of the coefficients 
of (1.1); i.e., 
1 71i(y) | 
) 1 T2(¥) |" 
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We have shown that 
F=U+iV and f= u‘\+ iv (2.3) 


are also 2-monogenic functions. We call them the 2-integral and the 2-deriva- 
tive of f respectively, and write 
- , af 
F= Sdxz, I va in Copia (2.4) 
zo dyz 
The definition of higher =-derivatives is obvious. We write 
fAiuaf, fraf, fraf', et. 
A comparison of the formulae (2.1)—(2.4) shows that 2-integration and 
>-differentiation are inverse processes. 
By applying 2-integration to constants we can obtain an unlimited num- 
ber of particular solutions of (1.1). The function f=1 is Y-monogenic, for 
u=1,v=0 isa solution of (1.1). Therefore 


° y dy 
Z“)(z) -{ Idsz = x+ if 
0 0 T1 


° - . ° Py , ° ° P 
is 2-monogenic, 1.e., “=x, v= J dy/t1 is a solution of (1.1). Next, 


es ae : y dy : v dy 
Z)(g) = 2 Zdsz = x? — 2 T2dy + 2ix — 
0 0 T1 0 T1 


is 2-monogenic and so are 


Z(z) = 3 f Z)dsz, Z% (2) = +f Z 4 d3z, - 
0 0 


It is not difficult to find a general formula for Z‘(z).1 We set Z =1 and 
define 


; t : v dy - y ; 
Y (y) = 1, Y()(y) -f Y(y) = 2f 72V “(y)dy, 
0 T1 0 


y 1 
Y®(y) = sf —YV(y)dy, --- 
0 T1 


Zs) = Do (% )arirvor(y), (2.6) 


v=0 


(2.5) 


Then 


1 Note that the superscript (m) does not denote differentiation. For the gas-dynamical 
equations written in the form (1.1), where 7; = 1, particular solutions corresponding to the formal 
powers Z) have been obtained independently of us by Professor S. Bergman (see footnote 14), 


and for the imaginary parts by Dr. A. Vaszonyi. 
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By repeated 2-integration of the function f=7, we obtain another set of 
~-monogenic functions, which we denote by 
t-Z(™)(z), 
We have 


i-Z(2) = id) (8) viemveo(y), (2.7) 
v=0 
where 
v . vy 
y*(y) = 1, V*?(y) -{ redy, Y* = 2f — Vt dy,--+. (2.8) 
0 0 T1 
Finally, fot any complex constant a=a+i8 we set 
a-Z™(z) = aZ(z) + BL i-Z™(z)}. (2.9) 
Clearly 
ds 
—~ ,a-Z(™) } = ga-Z'"—), (2.10) 


ayZ 


From the “formal powers” (2.9) we can construct new particular solu- 
tions. Obviously a “formal polynomial” of the mth degree 


f(z) = ado + ay-Z (cs) +--+ + a,-Z2™(s), a, ~ 0, 


is a L-monogenic function; i.e., its real and imaginary parts satisfy (1.1). 
It can be shown that there always exists a formal polynomial of the mth de- 
gree which takes prescribed values A;, Ao, +++, Anyi at m+1 prescribed 
points 21, > + *, Zn41. 

A “formal power series” 


is) 


f(z) = DY ay-Z™(z) (2.11) 
n=0 
represents a L-monogenic function provided the series converges uniformly 
and absolutely for sufficiently small values of | z|. It can be shown that any 
function which is =-monogenic around the origin can be represented in the form 
(2.11). The coefficients are given by a “Taylor formula”: 


nia, = f'")(0). (2.12) 


If in defining the formal powers we would extend the integration not from 
0 to z but from another fixed point 29, we would obtain another set of particu- 
lar solutions, which we denote by 


a-Z\)(z). 
bd) 
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Any 2-monogenic function defined in the neighborhood of 2, can be repre- 


sented in the form 


The same method can be applied to the more general equations 


O1(X)Ur = 71(¥)2%y, 
2.14) 
O2(X)Uy = — Tol Viz 


(o;>0,7;>0,7=1, 2). We denote the matrix of the coefficients of (2.14) by 2, 


4 
| 


and call a function f=u+iv =-monogenic if u and v satisfy (2.14). The 


2-integral, F= U+7V, and the =-derivative, f‘=u‘+iv‘, of f are defined by 


. 
* 4 u 
U = | ooudx — Tovdy, J -{ —dx + —dy, 
zg Fl 71 


e 


and 


U, V and wu‘, v‘ satisfy not (2.14) but the associated equations 


Uz 
= T] V)vy 
72 Yr) 
(2.45) 
Uy 
= — T2(V)V;z 
O1\ xX) 
Thus F and f‘ are 2’-monogenic, where 
1 
TI 
do 
a = 
1 
T2 
| 


However, 2”=2, so that the D’-integral of F and the &’-derivative of f* 


are 2-monogenic. 
It is clear in which way 
ers. We set 


we must change the definition of the formal pow- 
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Z” = f raze, ZZ = 2 f Z%des where Z( = f sass, 
Z3) = 3 [z dys where Z® = 2 f ze dxz,---, 


e 


and similarly for 7-Z“™. 

Remark. If in (1.1) 7, and rz are of opposite sign, the equations (1.2) are 
of hyperbolic type. The definition of =-monogenic functions, 2-integrals and 
derivatives, formal powers and formal power series remains the same. How- 
ever, it is mot true that « and v are necessarily analytic functions of x and y. 
Neither is it true that all =-monogenic functions can be represented in the 
form (2.13). 

The integrals defining Z‘" and i-Z‘ must not necessarily converge when 
T, OF T2 Vanishes at 2 or at z. If they do converge, they represent 2-monogenic 
functions. In this way it is possible to obtain particular solutions of partial 
differential equations which are of elliptic type in one part of the plane and 
of hyperbolic in another. 

3. Correspondence between >-monogenic and analytic functions. Let 


oo 


g(z) = PR G,2" (3.1) 


n=0 


be an analytic function of z. We define the =-monogenic function 


=a 


f(z) = Do an-Z™(2) (3.2) 


n=0 


and say that f corresponds to ¢ at the origin.’ 

It can be shown that series (3.2) converges in some neighborhood of the 
origin. If (3.1) converges everywhere, so does (3.2). If (at the origin) the 
~-monogenic functions f and g correspond to ¢g and y, then {+g corresponds 
to o+y, f'™ to go, and af to ag, a being a real constant. 

The concept of corresponding functions may be of use in discussing physi- 
cal problems. Assume that some physical phenomenon is described by equa- 
tions of the form (1.1). Often a simplifying assumption (for instance: the 
assumption of incompressibility) leads to Cauchy-Riemann equations. Sup- 
pose we possess an interesting or typical solution of the simplified problem. 
It will be given by an analytic function ¢=¢(z). We may expect that the solu- 
tion of the original more complicated problem given by 


¢ = f(z), f corresponding to ¢, 
will be of the same general character. 


? Thus the formal powers “correspond” to ordinary powers. 
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We define the 2-monogenic exponential and trigonometric functions (de- 


pending upon a real parameter a) 
E(a,2), S(a,z), C(a,2), i-E(a,z), t-S(a,2), 7i-C(a, 2) (3.3) 


as the functions which (at the origin) correspond to the analytic functions 
e*?, sinaz, cosaz, ie*%*, isinaz, 1COS az. 
We have, for instance, 


= en © (— 1)%q2t! 
E(a, z) = > Z‘™)(z), t-S(a, Zz) = pm : — 4-Z(2ntl)(z), 
n=0 i ! n=0 (2n + 1) ! 


o 


The function C(a, z)+i-S(a, z) corresponds to e‘***; E(a, z)+E(—a, 2), 
E(a, 2) —E(—a, 2) correspond to 2 cosh az and 2 sinh az. 
It is worthwhile to divide the functions (3.3) into their real and imaginary 


parts. A simple calculation shows that 


E(a, z) = e%*[c(a, y) + is(a, y)] 





S(a, z) = sin ax ch(a, y) + icosax sh(a, y) (3.4) 
C(a, z) = cosaxch(a, y) — isin ax sh(a, y) 
where 
2 =_ 1)*%a:2*! « q?zntl 
s(a, y) = >> rere YOtD(y), sh(a, y) = — — Y@at)(y) 
no (2+ 1)! nwo (2n + 1)! - 
(3.5) 
3 o.? P ° \ 
; 2 (— 1)*a* ia : en 
c(a, y) = - — Y(2")(y), ch(a, y) = —— Y(2")(y), 


n—o0 (2n + 1)! n-o (2n)! 


Note that the series (3.5) converge for all real values of y and for all complex 
values of a, so that we may write 
sh(a, y) = — ts(ta, y), ch(a, y) = c(ta, y). (3.6) 
Similarly we have 

i-E(a, z) = e%*[— s*(a, y) + ic*(a, y)] 

i-S(a,z) = — cosax sh*(a, y) + isin ax ch*(a, y) (39) 

i-C(a, z) = sin ax sh*(a, y) + icosax ch*(a, y), 
where the functions s*, c*, sh* and ch* are defined in the same way as in (3.5), 
Y™ being replaced by Y*™, 

The right hand sides of (3.4) and (3.7) have the form ,(x)¥i(y) 
+i,(x)W.(y). Thus we see that the real and the imaginary parts of the 
~-monogenic exponential and trigonometric functions coincide with the func- 
tions which we would obtain by solving equations (1.2) by the method of 
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the separation of variables. Therefore the functions (3.3) can be used in solving 
certain boundary value problems. 

The functions E, - - - ,t-C possess many properties analogous to those of 
ordinary exponential and trigonometric functions. For instance, the 2-differ- 


ential equations 


E. = aE, S* = al, C‘ = — aS. (3.8) 
hold. 
In the same way we could define correspondence at an arbitrary point 29 
and discuss the functions E(a, z), S(a,z), - - - , which at 2) correspond to the 
zo Zo 


analytic functions e*“-*, sin a(z— zo), °°. 

In some special cases the formal powers and the functions E£, S, C, can be 
expressed in a closed form by means of known functions. Such a case will be 
discussed in the next section. 

4. The case r;=72=y~. In this section we shall consider the special case 

| 1 yo? If 
> =| aa = (4.1) 


i] 1 y? ii 
For the sake of simplicity we assume that 
p20, 


and set 


p = 2% +1. (4.2) 


Along the real axis tr: =7T2 vanishes (except in the trivial case p=0). How- 
ever the integrals defining Y™ exist, and therefore so do the functions Z™, 
E, C, S. The integrals defining Y*™ converge only if <1. In this case only 
can we define the functions 7-Z™, «-E, i-C, and12-S. 

From (2.5) we obtain by a simple calculation 
, 


Y (+1)(y) —_ cE ee 


(4.3) 


| 
| 


Y()(y) = 


so that, by (4.3) and (2.6), 


> (- 1)” y 2v 
Z(s) = nIT(q + 1 — ! - gn (=) 
sii 4 ro vi(n — 2v)'T(qg +v + 1) ’ 2 


[(n—1)/2] _ 1 v y 2v+1 
oo iy? >» ——— ( _ ) voi(2) \. 
: vi(n — 2» — 1)! (¢g +r + 1) 2 


v=0 








Simpler expressions for Z™ will be found later. 








176 LIPMAN BERS AND ABE GELBART [Vol. I, No. 2 


Next we determine the 2-monogenic exponential and trigonometric func- 


tions. For the function s defined by (3.5) we get 


: (— 1)*1-3--- (2+ 1) 
sla, Vv) = >> - 7 qrtlyptir+l 
: (2y + 1)o+i1)--- (p+ 24+ 1) 
ae 1)°T gr 1) 
= z. a lyptar+l 
; Pile + » 4+ 2)2""* 
(— 1)T(q) 
= I "i l)bavy? > (a y)?, 
: y) : 


yIN(g +v +2 


Recalling the definition of a Bessel function* we see that 
S(a, y) = 297 qr l)a ay@tl] (ay). (4.4) 

Similarty we prove that 
Cla, y) = 297" gt l)a ty qa, ay). (4.5) 


For p=0 (i.e. g= —4) these formulae transform themselves into the familiar 
expressions for J;/2 and J_1/2. For in this case c(a, y) =cos ay, s(a, y) =sin ay. 


By (3.4), (3.6), (4.4) and (4.5) we have 


. - 1 { \ 
E(a, z) = 2°T(q + 1)a~%e* +4 JT (a) 7 ¢3 J ee aV) i, (4.6) 
. . - ° ° ° , ° ) wa 
S(a, 2) = 2°T (gq + 1)(ia)-*) y~? sin axJ, (tay) + ytt! cos ax q4i(tay)}, (4.7) 
~ ‘. ; ( ° . ° ) ’ 
C(a, z) = 2°T(g + 1)(ta)—*) y~? cos axJ,(iay) — y**! sin ax gii(tay)j. (4.8) 


In the case when 
p> 


the correspondence (at the origin) between analytic and 2-monogenic func- 
tions can be expressed by an integral formula. Since only the Z’s are de- 


fined we may consider only such analytic functions which are real on the real 


axis. 
For our case the equations (1.2) have the form 


(4.9) 
(4.10) 


yAu + pu, = 0, 
yh am pry = (). 


It is known that the solution of (4.9) which coincides with the analytic func- 
tion ¢(z) on the real axis is given by the generalized Laplace integral‘ 


3 See for instance, E. T. Whittaker and G. N. Watson, A course of modern analysts, 3d ed., 
Cambridge Univ. Press, 1920, p. 359. 
4 See H. Bateman, Partial differential equations of mathematical physics, Cambridge Uni- 


versity Press, 1932, p. 408. 
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1 r 
“u=- f g(x + ty cos a) sin?—'a da. 
nr Jo 
A simple argument shows that a solution of (4.10), connected with u by (1.1), 


and vanishing on the real axis is given by 


1 °. 
y= — 1yP —f g(x + iy cos a) sin’—!a cos a da. 
T 0 
Thus the =-monogenic function corresponding (at the origin) to ¢(z), is given 
by 
a i : , 
f(z) =- g(x + iy cos a)(1 + y” cos a) sin?—a da. (4.11) 
TJ 9 
In particular we have 
r 1 ” , . . 
Z™)(z) = - (x + iy cos a)"(1 + y? cos a) sin?—'a da. (4.12) 
T 0 


The integral (4.11) represents a 2-monogenic function even when ¢ pos- 
sesses a pole at z=0. We define 


1 r 
Z‘-*)(z) =- f (x + iy cos a)-"(1 + y? cosa) sin?'ada. (4.13) 
us 0 


It is easily seen that 


dy 
—ZO-"(z) = — nZO—"-1(z), 
dyz 
We now assume that p is an odd integer. In this case the formal powers 
can be expressed by Legendre polynomials. 
Introducing polar coordinates r, ¢, 
x =rcosg¢, y=rsing 


and recalling that the associated Legendre functions P? admit the representa- 


tions® 


1 2° (m+1) I'(n+m+1) . 


T T'(2m+1) r(n—m+1) 0 


n m m 


ry P,(cos ¢)= 





- n—m ., 2m 
(x+y cos a) sin ada 


we obtain from (4.12) 


“3 I'(p)T(n + 1) . i” 
Oe aa epee Lh te OP (4.14) 
2°1(q + 1)P(p + n) 


5 Bateman, |. c., p. 407. 
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Note that this formula is valid for an arbitrary p>1, if P?’, where m is not 
an integer, is defined by Hobson’s formula. ® 
If, as we assumed, g is an integer, we have by definition that 
g \ # (q) 
Prig(cos ¢) = sin g Pry,(cos ¢) 
where parenthesis on the superscript indicates differentiation, Then, by (4.13), 
(n) (p _ 1) !n! n_ (q) 


ReZ (zs) = —————r Pn,,(cos ¢). 
2%q!(2qg + n)! 


By (2.2) 
y? 0 
Im Z‘") = — aa ae Te grey 
n-+1 dy 
so that by (2.6) 
(n) (p — 1)!n! p n—1 (q+1) (q) 
Im Z ie y > A {COs y Pasqei(cos gy) — (1 + n) Pu4g+1(COs ¢)} 
21g i\p TT n)$ 
and 
(n), . (p—i1)ln! , (q) , 
Z (s)= , id i (p + n)Prs¢q(cos ¢) 
29q'(p +n)! 
oF q+1) (q) 
+ iy sin g[cos gPn4941(cos g) — (2 + 1) Pasq41(Cos y)}}. 


An analogous formula can be obtained for Z“”. 

If 2 has the form (4.1), particular solutions can be obtained also in an- 
other way. Introducing the polar coordinates r, g as new independent vari- 
ables, we obtain from (1.1) the system 


. ly, = SIN’ Vo 


r?—lyo = — sin~"¢ 2,. 
Denoting r+ig by ¢ and setting 


oe (4.15) 
25 : 4.15 


r?-!_ sin-?9 


we see that a 2-monogenic function of z is a Yx-monogenic function of ¢ and 
vice versa. 

Constructing Lx-monogenic formal powers, Z“(£), we obtain particular 
~-monogenic functions of z. 

For instance, if p=0 (Cauchy-Riemann equations), we have 


6 Whittaker and Watson, loc. cit., p. 325. 
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.-1, | 
~~ = 1] back 
ra 
and 
Zs = log z, i-Z, = i log z. (4.16) 


Of course, this method works not only in the case (4.1) but also when 
T, and 72 are of the form x*y'(x?+y?)™. 

5. Potential flow of an incompressible fluid with retational symmetry. In 
this and the following sections we shall briefly consider some mechanical ex- 
amples leading to equations of the form (1.1). 

Perhaps the simplest example is that of a rotationally symmetric potential 
flow of an incompressible perfect fluid.? We introduce cylindrical co-ordinates 
p, 0, z, p=0 being the axis of symmetry. Let qi, g2, g3 be the components of 
the velocity vector g in the direction of increasing p, 0, z. Then g2=0 and 
0q:/00 =0, dq;/00 =0. Using the well known formulae for div g and curl g we 
may write the continuity equation in the form 

i : @ 0g3 
div g = — — (pqi) + — = 0. 
p Op Oz 


The condition of irrotationality takes the form 


- 0q1 0q3 | 
| curl gj] = |— — —|=0 

| OZ Op | 
These equations imply the existence of a function u (velocity potential) and 
of a function v (Stokes’ stream function) such that 
gi = Up, q3 = U2, 


(5.1) 


ma te PI3 = Up- 


u-+iv may be called the complex potential of the flow. We have, by (5.1) 


1 
uz = —% 
p 
(5.2) 
1 
t = — —2D,. 
p 


In order to use our previous notations we shall write x, y instead of 2, p. 
Then we see that the complex potential u+iv is a =-monogenic function of 
x-+1y=2, where 

7 See, for instance, M. M. Munk, Fluid mechanics, Part I1, in W. F. Durand, Aerodynamic 
theory, Springer, 1934, vol. 1, pp. 260-263. 
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Since > is of the form (4.1) we may use the results of section 4. For the 
formal powers we obtain [cf. (4.2) and (4.14)], when g=0, x=r cos ¢, 
y=rsin ¢, 

/ ‘ii j . . , 
Z‘") = ——— j(n + 1)P,(cos ¢) + ty sin y|cos ¢ P.i41(cos ¢) 
n+ 1 
— T0 + 1) Pn+1(cos ¢) |}. 
Using a well known identity this may be written as 
r i 9 , ~ \ 
Z‘” = r"P,(cos ¢g) + r™-ly?P", (cos ¢). (5.3) 


n+ 1 


The real and the imaginary parts of the right hand side of (5.3) are the 
known polynomial solutions of the second order equations for u and v: 


yAu + u, = 0, 
0. 


(5.4) 


yAv — vy 
Now we introduce in the x, y-plane polar co-ordinates r, g (this amounts 
to introducing spherical co-ordinates in the physical space). Then u+vw will 
be a Ss-monogenic function of {=r+ig, where [cf. (4.15) | 
r? sino 
> — 


1 sin“ 9 


Forming the first 24-monogenic formal powers we get 


) ] 1 
5 1) */ ° - 
Ze =1- + 1(1—cosg) = — ( — + 400s °) + const. (5.5) 
. r r 
and 
—— 1 P - 
4-Z, = — log tan 3g +77. (5.6) 


ix/2 


In the case of a two-dimensional flow (Cauchy-Riemann equations) Z;’ was 
the complex potential of a source and i-Z{? that of a vortex [cf. (4.16) }. 
Similarly, in our case, (5.5) is recognized as the complex potential of a sink. 
(5.6) is the potential of a flow for which the streamlines are circles around the 
origin. The line y=0 is a singular line; it is covered by continuously distrib- 
uted sources and sinks. 

Other particular solutions can be obtained by forming higher 2»-mono- 
genic powers and by iterated 2-integration of (5.5) and (5.6). 





1943] DIFFERENTIAL EQUATIONS IN MECHANICS OF CONTINUA 181 


>-differentiating (5.5) we obtain (exactly as in the case of a plane flow) 
the complex potential of a doublet. In fact this 2-derivative equals 
cos ¢ _ sin*’g 


—i 
r? r 





Repeated 2-differentiation leads to doublets of higher order. It is easily seen 
that the 2-monogenic functions obtained in this way coincide with the formal 
negative power defined in section 4. In our case, (p=1), these powers have 
the form 


Zn = r~"—1P,. (cos ¢) —_ diya? "* PS (cos ¢). 


>-differentiation of 7-Z, does not lead to new particular solutions, since the 
>-derivative of (5.6) coincides with (5.5) (but for an additive constant). 

Finally we note the form of the 2-monogenic exponential and trigonomet- 
ric functions. From (4.2) and (4.6)—(4.8) we obtain, setting p=1, 


E(a, z) = e2*{ Jo(ay) + iyJi(ay)} 
S(a,z) = sin axJo(iay) + y cos axJ;(iay) (5.7) 
C(a, 2) = cos axJo(iay) — y sin axJ;(iay). 


6. Torsion of elastic bodies of revolution.* We consider an elastic body 
of revolution and introduce cylindrical coordinates p, 0, z, p =0 being the axis 
of symmetry. The physical components of the stress tensor (in the above 
coordinates) shall be denoted by a,x, those of the displacement vector by u,, 
i, k=1, 2, 3. Because of the symmetry these quantities do not depend upon @. 
Furthermore, we assume that 


u; = us = 0. (6.1) 


Then all o;, vanish except 012 and o23. The condition of equilibrium takes the 


form 


0012 0023 2o12 
ieee i samen iis same iy A (6.2) 
Op 02 p 


whereas the stress-strain connection is given by 


0 {Ue Ou2 
ous -c—(*), oe: (6.3) 
Op\p Oz 


8 See A. E. H. Love, A treatise on the mathematical theory of elasticity, 4th ed., Cambridge 
University Press, 1934, pp. 330-332, and A. Timpe, Die Torsion von Umdrehungskorpern, 
Math. Annalen 71, 480-509 (1912). 
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Substituting from (6.3) into (6.2) we obtain 


07u3 O72 1 Ou2 u 
- —— — = @, (6.4) 


Op? Oz” p p p* 


If the function y is defined by 
y = pu, 


(6.4) may be written as 
1 
Yoo t¥ze2 —-—¥, = 0 
p 


This equation implies the existence of a function g connected with by the 


equations 
1 
,=—y 
‘i p 3 (6.5) 
1 
= ——v 
p 


yg and y can be interpreted as the velocity potential and Stokes’ stream func- 
tion of a rotationally symmetric potential flow [cf. (5.2)]. The particular 
solutions of (6.5) have been discussed in the preceeding section. 

Let our body be a cylinder subjected to a deformation of the kind con- 
sidered above. If the displacements on the boundary are given, the displace- 
ments in the interior can be found by integrating (6.4) under the boundary 
condition: 

¥=x(z), for p= + po, 
po being a constant and x a given function. Plainly, this boundary value 
problem can be solved by representing x by a Fourier series (or by a Fourier 
integral) and by using the particular solutions (5.7). Similarly the solutions 
(5.7) can be used if the boundary values are prescribed on the lines z=const. 
In this case developments in series of Bessel functions should be used. The 
discussion of convergence will be found in the quoted paper of A. Timpe. 

Equation (6.2) can also be written in the form 

0 0 
(p?o12) + — (p*o23) = 0. 


Op Oz 


This equation can be satisfied by introducing a stress function v and setting 


2003 = — (6.6) 
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On the other hand, writing 





Guz 
“= ; 
. - p 
we obtain from (6.3) 
Ou Ou = 
te 94 7 = > * (6.7) 
Op Oz 
Comparing (6.6) and (6.7) we get 
1 
u, = 3 Vp 
(6.8) 
1 
up = — — Pz. 
p? 


I'he torsion problem consists of integrating equations (6.8) under the boundary 
condition 


v = const. on the boundary. 


his condition expresses the fact that the surface of revolution is free of 
stresses. 

In order to maintain our previous notations we shall write x, y instead 
of z, p. Then (6.8) shows that u+7v is a -monogenic function of x+7zy =z, 
where 


I 
3s . | | 
~ 
! 
i 


_ 
e 
wo 


lhe second order equations for u and v take the form 
yAu + 3uy, = 0, (6.9) 
yAv — 30, = 0. (6.10) 
Again » has the particular form discussed in section 4. From (4.2) and 
(4.14) we obtain, setting p=3, 


Nn. 


>) , 
Z‘*)(z) = r™) (n + 3)Prsi(cos ¢) 
3)! 


v7? 


° « . , 
+ iy? sin y|cosy P,42(cos ¢) — (nw + 1) Pry2bcos ¢)]} 
(r, g are polar coordinates). This is easily transformed into 
1 


Z™) = et ; 2) {r"P...1(cos ¢) + irnt? sinty P.41(cos y)}, (6.11) 
(m + 1j}(n y 
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The imaginary parts of (6.11) are the known polynomial solutions of (6.10).° 
Next we form the 2-monogenic exponential and trigonometric functions. 


From (4.6)—(4.8) we obtain 


2 1 
E(a,z) = — cor} —Ju(ay) + iy'Juay)} 


a y 
2(1 

S(a,z) = —_~ sin axJ(iay) + y? cos a(iay)\ 
ay 
2(1 

C(a, z) = —4{— cos axJi(iay) — y* sin aaJa(iay)\ ; 
aly 


The imaginary part of E again coincides with a well known particular solution 
of (6.10). '° 

If we introduce the polar coordinates as new independent variables, u+-iv 
becomes a 2%-monogenic function of f =r+7g, where 
sean r* sin~* ¢ || 


“~~ — 


I} r? sin-* ¢ || 


Forming the first 2%-monogenic formal powers we obtain 


1) 1 ° 
Zz. =1 — 3 + if — 3 cos ¢ (sin*e + 2)] (6.12) 
r 
ight cos 9 © as 
‘4 = ———— — , log tan 39 + t—. (6.13) 
in/2 2 sin? g 3 


(6.12) describes the torsion of a cone with the vertex at the origin. (6.13) 
represents the torsion of a sphere around the origin, the torque being trans- 
mitted through a singular line along a diameter. 
2-differentiating (6.12) we get 
dz a) cos © _ sinty 


ct te -— 9 
d yz r’ Tr 


(6.14) 








The equations (6.8) can also be interpreted hydrodynamically. u—iv can 
be understood to be a complex potential [(velocity potential) +i(Stokes’ 
stream function) | of a rotationally symmetric potential flow in a five dimen- 
sional Euclidean space. (This interpretation is due to Arndt).!! If we adopt 
this point of view, Z again represents a source, and its 2-derivative, (6.14), 


* Love, I. c., p. 331. 


10 Love, l.c., p. 332. 
11F,C. Arndt, Die Torsion von Wellen mit achsensymmetrischen Bohrungen und Hohlréu- 


men, Thesis, Géttingen, 1916. 














a erat ate 
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a doublet. Combining “flows” of the form (6.12) and (6.14) Arndt solved the 
torsion problem for a number of bodies with holes. 
Repeated 2-differentiation of (6.14) leads to doublets of higher order. 
These complex potentials are equivalent to the Z‘”’s defined in section 4. 
7. Two dimensional potential gas flow. Let p and p represent the pressure 
and density of a perfect compressible fluid.’? If heat conduction is neglected 
the above quantities are connected by the relation 


p= (po/po)p 


where y is the ratio of the specific heat for constant pressure to the specific 
heat for constant volume, the subscript zero referring to the fluid at rest. 
Introducing the local velocity of sound, a, given by 


we have by Bernoulli's equation 
2 


a = ao — [(y — 1)/2]q, 


where g is the magnitude of the velocity, and 


Y- 1 g? 1/(y-1) 
p = »s[ 1 — <| e 
2 aR 


The continuity equation for a steady flow has the form 


div (2 i) = 0, (7.1) 
Po 


where g is the velocity vector. If the flow is irrotational, 


curl’g = 0. (7.2) 


For a two-dimensional flow, g = (q cos 8, q sin 8, 0), @ giving the direction of the 
flow. (7.1) and (7.2) then imply the existence of two functions, ®(x, y), 
W(x, y), such that 


gcos0= 9, qsin@ = ,, 
and 
Po . Po 
gcos# = —YW,, gsin@= — —W,. 
p p 


12 See, for instance, G. I. Taylor and J. W. MacColl, The mechanics of compressible fluids, 
in W. F. Durand, Aerodynamic theory, Springer, 1935, vol. 3, pp. 229-230. 
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Thus the velocity potential ® and the stream function WV are connected by 
the equations 


Po 
’,=—YW, 


Po 
b= ——W,. 
p 
Since p is a non-linear function of the derivatives of the unknown functions 
® and YW, the above equations are non-linear. However, linear equations can 
be obtained by a transformation due to Molenbroek" and to Chaplygin;" 
namely, by introducing as new independent variables the quantities @ and gq. 


The equations then become 


9 


Po 
= — qv, 
P (7.3) 

Po q* 1 
6,= ——(1-—)—W. 

p a*/} q 

To simplify the formulae Chaplygin introduces a new variable 
2(y — 1) 


se 

In what follows we write x, y instead of 0, ¢. (x, y are not, as they were before, 
the coordinates in the physical plane). Then (7.3) shows that the complex 
potential +7 is a 2-monogenic function of x+7y, with 


2y | 
a ‘ai y)™ | 1 
‘- 1 — (2m + I)y 


2y(1 a y) l+m | 


i} 


We see that 
7, = 1/72 =0 for y=0O, t2=0 for y=1/(2m+ 1). (7.4) 


13 P. Molenbroek, Uber einige Bewegungen eines Gases mit Annahme eines Geschwindig- 
keitspotentials, Arch. d. Mathem. u. Phys. (2). 9, 157 (1890). 

4 A. Chaplygin, On gas jets, Scientific Annals of the Moscow University, Section for Math. 
and Phys, 21, 1-121 (1904) (Russian). Cf. Th. von Karman, Compressibility effects in aerody- 
namics, Journ. Aeronaut. Sc. 8, 337-356 (1941). (This paper also contains an extensive bibli- 
ography.) For a detailed description of the Chaplygin transformation in English see 
S. Bergman, The Hodograph Method in the Theory of Compressible Fluids (suppl. to K. O. 
Friedrichs, R. v. Mises, Fluid Dynamics), mimeographed lecture notes, Brown University, 1942. 
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The formal powers a-Z,,“ cannot be defined for Im z9=0. For 
0 < y < 1/(2m + 1), (7.5) 


i.e. for subsonic flow, we have the elliptic case, 7; >0, 7 =1, 2. 

We postpone a systematic discussion of the application of our method to 
this case and shall show only how an important set of particular solutions 
(due to Chaplygin) can be represented in terms of our functions. 

The second order equation for Y takes the form 


1 — (2m + 1)y 0 
; ai - ee V..+ a {2y(1 — y)-™¥,} = 0. (7.6) 
Chaplygin gives as particular solutions of (7.6)% 
¥ = yFq(y) cos ax, 
x = yF.(y) sin ax, 
where a>0O and F,(y) is the hypergeometric function 
F.(y) = F(aa, ba, 2a + 1, y) 
the constants da, 0. being determined by the conditions 
da + b, = 2a — m, Gaba = — ma(2a + 1). 
Consider the =-monogenic function 
H=¢+ wy; 
we have 
H = 9.+ Ww. = 2y1— y+ Wz 


2y(1 — y)-™yt [ak a(y) + yF d (y)] cos ax — iay*F(y) sin ax. 


I 


Therefore 


h= H', 
satisfies the differential equation 
7h P 
—-+a*h = 0, 
Ox? 
i.e. the 2-differential equation 
h + a®h = 0. (7.7) 
8 Chaplygin, 1. c. Similar solutions have been considered by F. Ringleb, Exakte 


Lisungen der Differentialgleichungen einer adiabatischen Gasstrimung, Zeitschr. angew. Math. 
Mech. 20, 185-198 (1940). The method by which Chaplygin solved the jet probiem and Ringleb 
computed the complex potentials of a compressible doublet is actually identical with the use of 
the “correspondence” discussed in section 3. 
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Let yo be a fixed real number, satisfying (7.5). Then 


h(iyo) = aA, h‘(iyo) = ia®B (7.8) 
where 
1 m a—! : - 
A = — 2yo(1 — yo) yo [aFa(yo) + yola (yo) ] 4 
Qa (7 F 9) 
B= — yo Fa(yo)- 


It is seen easily that the solution of (7.7) with the initial conditions (7.9) is 
AC (a, z) + Bi-S (a, 2). 
tivo tvo 
This function is therefore identical with h. 2-integrating and suppressing a 
non-essential constant we obtain 


H =AS (a, z) — Bi-C (a, 2) 


Yo tVo 
so that 
y = In[A S (a, z) — Bi-C (a, 2z)]. 
ivo ivo 
In a similar way we obtain 
x = — Im[AC (a, z) + Bi-S (a, 2)]. 
ivo ivo 


It is interesting to note that if a is a positive integer, H=g+iy and G=w+ix 
are regular at the origin and that all their =-derivatives vanish at this point. 
This is possible because of (7.4). 
We also note that the first formal powers 
Z" and 4-Z°) 


represent the complex potentials of a compressible source and vortex respec- 


tively. 





16 Cf, Ringleb, |. c. 
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BOOK REVIEWS 


Mathematics of modern engineering. Volume II (Mathematical Engineering). 
By Ernest G. Keller. (One of a series written in the interest of the Ad- 
vanced Course in Engineering of the General Electric Company.) John 
Wiley and Sons, Inc., New York, 1942. xii+309 pp. $4.00. 


This book is the second in a series of three, the first of which was written 
jointly by the present author and R. E. Doherty. Roughly speaking, the pres- 
ent volume treats problems leading to ordinary differential equations while 
the third volume, apparently, will be devoted to applications of partial differ- 
ential equations. The material of the book is said to have evolved out of re- 
search engineering work and the Advanced Course in Engineering of the 
General Electric Company. It was prepared for that course and for graduate 
engineering work in general. It was decidedly not intended either as a text- 
book in mathematics or in engineering. Its main purpose is to illustrate, by 
means of many typical examples, the methods of what the author describes 
as a new phase in engineering and calls Mathematical Engineering. 

As in Mathematical Physics, these methods are essentially mathematical 
in character. Still, mathematical engineering is far from being a part of mathe- 
matics. Indeed, solving a practical problem consists of two parts. The formu- 
lation of the problem and its reduction to a form in which it lends itself to 
mathematical treatment constitutes by far the most difficult part. It requires 
not only knowledge, but also originality and instinct. Once the problem is 
reduced to a canonical form, its solution may be difficult but can usually be 
achieved by a well established technique. It should be a comfort to the non- 
mathematical engineer and to the mathematician alike that, with the numer- 
ous engineering examples of the book under review, the second step in each 
case reduces to routine. In presenting mathematical theories and methods the 
author has not felt tempted to deviate the least from common usage, thus 
indirectly paying tribute to the usefulness of “pure” theory. 

The great interest, and value, of the book lies in its large collection of ex- 
amples, taken from the rich experience of the author, and illustrating the 
application of mathematics to engineering problems. Some of the problems 
are given as simple exercises (without solutions), others as more or less simple 
illustrative examples to the mathematical text. The most valuable material is 
represented by examples where the reader can follow the whole procedure of 
formulating, and solving, difficult problems of engineering practice. For ex- 
ample, with an intricate problem concerning electric locomotives the factual 
information is presented in the form in which it once actually confronted the 
author. The setting up of the problem is done with masterly skill and the final 
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formulation leads to a surprisingly simple mathematical question. This is an 
example of mathematical engineering at its best. 

The style of the book seems to indicate that the author himself was inter- 
ested only in presenting this new and interesting material. As a matter of 
fact, the general theory is treated in a rather unsystematic way. The author 
carelessly mentions here and there what the needs of the moment require, 
sometimes altogether forgetting the reader. Apparently, the latter is not even 
expected to be familiar with vectors and matrices. On the other hand, the 
short summaries of much more advanced topics cannot seriously be meant as 
a first introduction. Occasionally the intended elementary character of the 
book is given up, and once (p. 197) the author almost apologizes for mention- 
ing problems which have already been answered by engineers. One example 
of the haphazard way in which theory is presented may suffice. Ordinary 
linear differential equations are first presented as they unfortunately used to 
be presented in elementary textbooks, without matrices. In connection with 
G. Kron’s theories more than twenty pages are devoted to matrices, and here 
the said equations are put into matrix form (with the danger of producing the 
false idea that this is an application of the somewhat mystically formulated 
first “generalization postulate of G. Kron”). Finally, in the last chapter, for 
reasons unclear to the reviewer, the same linear equations are treated among 
non-linear problems, this time in the more efficient matrix form. Now the 
author also gives good advice on how to compute characteristic values while 
during the first exposition he suggested impracticable methods. More trivial 
examples of the same nature appear frequently. Thus, treating the two- 
dimensional automobile (p. 21), the author says: “Evidently, by Koenig's 
theorem ....” For the theorem itself the reader is referred to p. 44, where 
he finds only a formula and no theorem. For, only when the author needed 
the theorem for the second time, did he formulate it in the course of the dis- 
cussion of a special problem. 

The book consists of three parts. Chapter I, “Engineering Dynamics and 
Mechanical Vibrations” treats applications of classical dynamics. The subject 
of vibrations seems most congenial to the author. Accordingly, special ex- 
amples are most numerous in the first chapter; and, on the whole, this chapter 
seems the best. Chapter II is devoted to an “Introduction to Tensor Analysis 
of Stationary Neworks and Rotating Electrical Machinery.” This is actually 
an introduction to G. Kron’s theories. In the reviewer's opinion this chapter 
is by far the most difficult one and certainly no introduction in the usual 
meaning of the word. The exposition is less personal, the examples less numer- 
ous and less illustrating. Chapter III is entitled “Non-Linearity in Engineer- 
ing.” It treats briefly of various methods of obtaining solutions, or approxi- 
mate solutions, of ordinary differential equations. Moreover, elliptic and 
hyperelliptic functions and integrals are treated and examples of their ap- 
plication to vibration problems are given. There is an extensive bibliography. 
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but altogether too many of the books mentioned will prove valueless (or 
confusing) to the average reader. 

These deficiencies are not serious enough to endanger the usefulness of 
the book. They are mentioned in the hope that their removal in a future edi- 
tion will render the book still more enjoyable and still more useful as a text- 
book. Even in its present shape the book should convince many a sceptical 
engineer of the power of mathematical reasoning when applied by a man who 
combines mathematical knowledge with a good engineering mind. 

W. FELLER 


Microwave transmission. (International Series in Physics.) By J. C. Slater. 
McGraw-Hill Book Company, Inc., New York and London, 1942. x+309 
pp. $3.50. 


This new book in the “International Series in Physics” deals very com- 
pletely with the theory of transmission of microwaves (electromagnetic waves 
with wave length of the order of magnitude of 1 cm.). Problems concerning 
generation and reception of microwaves are omitted except for a few brief 
remarks, because, as the author points out in the introduction, very little is 
known from a theoretical point of view about their actual operation. 

The treatment of the subject is based largely on electromagnetic theory, 
although the first chapter is devoted to conventional transmission lines for 
radio waves. The ideas developed there are kept in mind and used for pur- 
poses of comparison in the later chapters. Considerable use is also made of 
analogies with acoustics and optics. The former are particularly pertinent as 
the wave length of acoustic waves is of the same order of magnitude as that 
of microwaves. Much of the terminology introduced is taken by direct anal- 
ogy from the corresponding acoustic case. 

The chapter on the fundamentals of electromagnetic theory is fairly ele- 
mentary but uses vector notation. Meter, Kilogram, Second (Georgi) units 
are used throughout. A complete discussion of electrodynamics would of 
course be out of place in a book of this type. This theory is developed in suffi- 
cient detail, however, so that the subsequent chapters should be intelligible 
to any senior in physics or engineering. 

After the preliminaries on ordinary transmission lines and electrodynam- 
ics, the author goes on to consider in detail the theory of various types of 
transmission lines. Wave guides and coaxial cables are thoroughly discussed. 
Particular emphasis is placed on phenomena which occur when the cross 
section of a wave guide of cable is varied and which take place at junctions 
of lines of different characteristic impedance. Problems of matching imped- 
ances and the applications are included. The various types of waves that may 
arise are classified and their propagation characteristics are studied. The 
theory of radiation from antennas is given in detail, with an entire chapter 
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devoted to directive devices for antennas. The final chapter deals with the 
coupling of coaxial lines and wave guides. 

The author’s treatment of these subjects is extremely good, and his pres- 
entation is clear and logical with all extraneous matter eliminated. The book 
should find considerable use both as a text for the student and as a reference 
work for the specialist. It represents, moreover, a real contribution, in that 
its sets forth in concise and unified form the recent advances in microwaves— 
an extremely important field at the present time. It is an excellent book, to 
be recommended to all the scientists interested in this field. Those who wish 
more detailed theoretical calculations should be reminded of Stratton’s Elec- 
tromagnetic theory (McGraw Hill, 1941) and E. V. Condon’s Principles of 
microwave radio (Reviews of Modern Physics, 14, 341-389 (1942)) where a 
number of very interesting examples are discussed in full detail in view of 


practical applications. 
L. BRILLOUIN 


Introduction to non-linear mechanics. By N. Kryloff and N. Bogoliuboff. 
A free translation by Solomon Lefschetz of excerpts from two Russian 
monographs. (Annals of Mathematics Studies, No. 11.) Princeton Uni- 
versity Press, Princeton, 1943. iv+105 pp. $1.65. 


The major part of this monograph is a very condensed version of a Russian 
book by the two authors. By making a part of the work of Kryloff and 
Bogoliuboff accessible to American scientists, the monograph accomplishes 
an extremely useful purpose and will, no doubt, be welcomed by any one inter- 
ested in the subject of Non-Linear Mechanics. 

The monograph deals mainly with the equation 


d*x dx 
— +s = o(« =) (1) 
dt 

(for small values of €) and with the problem of finding approximate solutions 
of (1), free from secular terms and satisfying (1) to within terms of order e” 
with prescribed m. Problems studied in the monograph refer, generally speak- 
ing, to non-conservative physical systems and differ in this respect from the 
corresponding problems in the astronomical theory of perturbations. 

The mathematical tools are developed mainly from the viewpoint of ap- 
plications in Physics and Engineering, and mathematical rigor is not the pri- 
mary consideration. It should also be observed that the qualitative geometrical 
methods of Poincaré-Bendixon are not discussed in the monograph which 
deals exclusively with analytical approximation methods. 

The first chapter contains various examples in Mechanics and Electrical 
Circuit Theory which lead to non-linear equations of type (1). The theory of 
the first approximation is studied in Chapters II and III. Higher approxima- 
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tions are discussed in Chapter IV. The method of “linearization” described 
in Chapters V and VI consists in substituting for equation (1) a linear equa- 
tion which, in the first approximation, has the same solution as (1). The 
important thing is that this linear equation can usually be set up directly 
from the physical data of the problem, without reference to the exact equa- 
tion (1). Chapter VII discusses multiply periodic systems. Periodic disturb- 
ances are studied in Chapter VIII. This leads to equations of the type 


d*x dx 
m—-+ kx = 7G Z, =) 
dt? dt 


where the function f is almost periodic with respect to the variable ¢. The final 
chapter contains several excerpts from a monograph by the two authors ex- 
hibiting relations between the properties of the exact solutions and those of 
the approximate solutions. 

Throughout the monograph mathematical theories are illustrated by nu- 


merous examples of physical systems. 
W. HuREwIcz 


Heat transmission. By William H. McAdams. (Sponsored by the Committee 
on Heat Transmission, National Research Council.) Second Edition. 
McGraw-Hill Book Company, Inc., New York and London, 1942. 
XIV+459 pp. $4.50. 

The knowledge of the processes of the heating and cooling of bodies has 
gradually accumulated from the dawn of history as a mass of empirical facts. 
Since the time of Newton who expressed some elementary ideas on the laws 
of cooling, the thermal aspects of the world have been divided into the 
equilibrium aspects (thermodynamics) and the non-equilibrium aspects (heat 
transfer). 

All the phenomena of thermodynamics were found to be completely cov- 
ered by a few sweeping general laws which serve to correlate precisely seem- 
ingly unrelated phenomena. No corresponding simple generalizations have 
been found for the non-equilibrium thermal phenomena. The basic physical 
laws of conduction and radiation were discovered and applied during the last 
century and it was clearly understood then that the great majority of heat 
transfer problems encountered in practice could be adequately treated by 
these laws only if the attendant convective movements of a fluid were prop- 
erly accounted for simultaneously. It is the presence of these convective 
movements which has and still does retard the analytical development of 
heat transfer knowledge as compared with other aspects of science. 

The first edition of Heat Transmission by Wm. H. McAdams, Prof. of 
Chemical Engineering, Mass. Inst. of Tech., appeared in 1933. The emphasis 
of the book was on the presentation, for ready reference, of all the best 
empirical data and computational methods required for the solution of prac- 








194 BOOK REVIEWS [Vol. I, No. 2 


tical heat transfer problems. That the time was ripe for the general collection 
and presentation of these heat transfer facts is indicated by the fact that the 
book was sponsored by the Committee on Heat Transfer of the National 
Research Council. That the book was successful as a text and reference work 
is attested by the fact that in the 10 years between the first and second edi- 
tions 9 impressions were needed to satisfy the demand and approximately 
eight thousand copies were sold. 

One quarter of the first edition was devoted to conduction and radiation 
heat transfer, the physical laws, the physical constants of materials, the 
empirical and mathematical methods for the solution of problems. The re- 
maining 3/4 of the book presented information needed for the solution of 
problems involving convection. A chapter on Dimensional Analysis presented 
the basis of the correlation of empirical data. A chapter on Flow of Fluids 
presented laminar and turbulent friction and pressure loss data. Succeeding 
chapters dealt in turn with, the concept of thermal resistance, convection 
heat transfer by fluids inside pipes, fluids outside pipes, condensing vapors, 
and boiling liquids. The book closed with some 75 pages of Appendix dealing 
with thermal properties of materials, Bibliography and Subject Index. 

Like most branches of science, the field of heat transfer is one in which 
developments are being made at an ever more rapid pace. Ten years of work 
made a considerable revision of the first edition necessary if the book was to 
maintain its position of excellence in the field of heat transfer literature. 
During these years, as before the first edition, Prof. McAdams followed care- 
fully the developments as they were made and indeed was himself responsible 
for some of them. 

The new edition has again been sponsored by the Committee on Heat 
Transfer of the National Research Council. 

The book is written in the same spirit as the first edition and follows the 
same general outline described above. The subject matter of every chapter 
has however been completely revised. The thoroughness with which the new 
work in heat transfer has been reported is well indicated by the fact that the 
former Bibliography of 414 references has been expanded to 7839, i.e. almost 
doubled. These are not useless references either, as each one is followed by a 
page number of the book where it is mentioned or where its subject matter is 
pertinent. The new data has been so well assimulated into the old that the 
book has been increased by only 100 pages. 

The opening chapter of the new edition is a four page introduction in 
which the definitions of conduction, radiation and convection are given to- 
gether with the basic equation for each. A most important addition is made 
on page 4 where a table of “Approximate Range of Values of Heat Transfer 
Coefficients Ordinarily Encountered” is given. 

Heat transfer by conduction is covered in 20 pages and as in the earlier 
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edition presents only the most directly useful facts, in good part through the 
ideas of a mean area or a thermal resistance. Analytical methods of solution 
are hardly more than mentioned in passing although some analytical results 
in the form of resistance formulas are tabulated. 

The 17 page chapter on the Heating and Cooling of Solids remains much 
the same as before. Analytical results, a single sketchy derivation of one case 
is given, are presented as dimensionless charts, several new ones being in- 
cluded. A graphical method of solving transient heat conduction problems 
has been added and illustrated by an example. 

Radiant Heat Transmission has again been written by H. C. Hottel, a 
colleague of McAdams on the Chemical Engineering Staff of Mass. Inst. of 
Tech. The subject matter, 40 pages of it, has been completely revised. The 
matter of spectral distribution of energy and its effect on heat transfer, hardly 
more than mentioned in the first edition, is treated in as much detail as the 
accumulated data on absorptivity and emissivity of bodies will permit. Fol- 
lowing the derivation of the general differential equation and a sample inte- 
gration is a discussion of geometrical and “non-black body” factors. Con- 
siderable new results have been added here in the form of what might have 
been called theorems regarding interrelations between the factors under vari- 
ous circumstances. The radiant heat transfer from gases has likewise been 
completely rewritten. 

The chapter on Dimensional Analysis in the first edition was easily the 
weakest in the whole book. Considering the fundamental position which this 
type of analysis holds in the correlation of all convection data, it is important 
that the fundamental ideas as well as its method of use be made very clear 
to the reader. The new writing of this chapter is enormously better, in spite 
of what the reviewer regards as incorrect statements in connection with the 
“dimensional constants.” It is indicated that dimensional constants must be 
added if more than four dimensions (say M, L, 0, T) are used. Why four? 
Why not three (7 = ML*°0-? by kinetic theory)? The answer is simple, kinetic 
theory as such does not play any part in the heat transfer correlations con- 
templated in the book. For a similar reason and with the same advantage it 
is sometimes proper and desirable to use five or six dimensions without adding 
dimensional constants. 

The discussion of Flow of Fluids remains much the same as before except 
for the addition of some of the dimensionless correlations of velocity distri- 
bution in turbulent flow in pipes. The sameness, by chance, even extends to 
the page numbers for Chapter V; first edition pages 99-134, second edition 
pages 99-132. 

Once the idea of thermal resistance (or conductance) is introduced to 
permit calculation of series and parallel thermal circuits, the most difficult 
problem of analysis of heat exchanger performance is the question of “mean 
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temperature difference.” The latest pertinent facts are presented in Chapter 
VI together with a brief discussion of the important experimental problem 
of the measurement of surface temperature. 

The next four chapters deal with the details of the correlated heat transfer 
data on fluids inside tubes, outside tubes, condensation and boiling. The re- 
visions and additions are extensive and too numerous to mention in detail. 
They involve new data used to improve the accuracy of old correlations, 
new improved methods of correlating old and new data. New theoretical re- 
sults are mentioned where of practical importance as, for example, in the 
analogy between heat transfer and friction. Even some new types of heat 
transfer problems have found an important place, as for example film boiling. 

A new last chapter has been added entitled “Applications to Design.” In 
this chapter, structural details, economic considerations and the choice of 
“optimum” operating conditions are discussed. The criterion of “optimum” 
is the fixed charges plus operating costs which depend upon the power re- 
quired to pass the fluids through the heat exchanger. Differentiation gives 
the conditions for minimum cost which are solved for the optimum conditions 
sought. 

Throughout the book a generous number of illustrative examples are 
given. As a further aid to the reader, each chapter is concluded with a page 
of problems on which to test one’s real understanding of the subject matter. 
The new edition of Heat Transmission assures the continuation of the 


high reputation of the earlier work for another decade or so. 
HOWARD W. Emmons 


Empirical equations and nomography. By Dale S. Davis. McGraw-Hill Book 
Company, Inc., New York and London, 1943. ix+200. $2.50. 


The first part of the present book will prove useful to engineers and sci- 
entists confronted with the task of finding “a simple mathematical expression 
fitting some experimental data.” The first two chapters are devoted to a dis- 
cussion of frequently used functions which contain a certain number of con- 
stants. Methods of determining the numerical values of these constants are 
given. The third chapter is concerned with the analysis of experimental data 
involving two independent variables. It is obvious that this problem cannot 
be treated comprehensively in the space of some twenty pages allotted to it; 
the reviewer believes, however, that the very useful graphical method de- 
veloped by A. Lafay [Génie Civil, 40, 298 (1902) ] should have been discussed. 

The second part contains an introduction to Nomography. Chapters IV 
to IX bring the most important types of alignment charts for three and more 
variables and Chapter X is devoted to special slide rules. 

Throughout the book the mathematical argument is of an elementary 
character; the discussion is illustrated by numerous examples, a large portion 


of which are concerned with problems of chemical engineering. 
W. PRAGER 
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Electromagnetic Waves. By S. A. Schelkunoff. D. van Nostrand Inc., New 
York, 1943. xv +530. $7.50. 


This new addition to a well-known series has been awaited with much in- 
terest by all those acquainted with Dr. Schelkunoff’s contributions to propa- 
gation theory, and it will be found that.their expectations have been entirely 
fulfilled. This monumental piece of work is equally remarkable for the origi- 
nality and consistency of its approach as for the wealth of information con- 
tained in its five hundred densely packed pages. 

The author's systematic use of the harmonic oscillation, with complex 
variables and coefficients, is in line with the marvelous development which 
has occurred in the communication field during the last fifty years. Alternat- 
ing current theory, then acoustics, then vibrational mechanics successively 
dropped the differential equations which physics offered as a. basis and sys- 
tematically restricted themselves to harmonic oscillations. This has resulted 
in the replacement of the differential operator by iw, leading to a tremendous 
simplification of steady-state analysis, which has been reduced to the calcula- 
tion of amplitude ratios and phase differences. The genuinely difficult prob- 
lems have not disappeared for all that but are now relegated to Fourier or 
Laplace transform theory, and it has become apparent that an enormous 
field of application can be covered by purely algebraic processes. 

Not the least advantage of this method has been the unification brought 
into the three chapters of technical science mentioned above. Electrical 
impedances gave the model after which acoustical and mechanical impedances 
were fashioned; and mixed mutual impedances, thereafter, made it possible to 
write the equations of electro-mechanical or acoustico-mechanical transducers. 
There was an exciting era of intense development in this field during the 
twenties; and it was amusing to hear at that time, and even a good deal later, 
irate die-hards denouncing “impedances” with bitter irony or viewing with 
alarm the spread of “analogies.” 

Dr. Schelkunoff has set about to carry this point of view into Electro- 
magnetic Theory, and it may well be that his will be the honor of having 
brought into the fold of harmonic oscillation theory the last chapter of 
Physics which still had to be incorporated. (One might think of Optics, but 
of course half of the book is really Optics.) Having given, in the first pages of 
his fourth chapter, a short and quite personal derivation of Maxwell’s equa- 
tions (1-15, p. 69), Dr. Schelkunoff without taking breath adds immediately: 
“Since we are concerned primarily with fields varying harmonically with time, 
we replace the instantaneous field intensities and current densities by the 
corresponding complex variables and write Maxwell’s equations as follows: 


fenas = — fiouttas ~ f M,dS, 
fataas f fi (g + iwe)EndS + f ras. 


(1—16) 








198 BOOK REVIEWS [Vol. I, No. 2 


Thus the sacrosanct Maxwell equations are swept away with movie-like 
swiftness, and instead we have the steady-state equations of a medium char- 
acterized by a distributed series impedance twy and a distributed shunt ad- 
mittance g+iwe (p. 81). 

The analogy with a transmission line whose series inductance is wu, shunt 
conductance g and shunt capacitance e¢, all taken per unit length, is inescap- 
able (p. 243). In particular the above primary constants simply beg to be 
transformed into the familiar secondary constants of transmission line theory; 
here the intrinsic propagation constant o and the intrinsic impedance 7 are 


defined by 


/ . 
/ Wop 


o = Vieng + iwe), n="/ sa Ny (9-1) 
WE 


(p. 81) (o is in neper/meter, 7 in ohms; the book is written in MKS—p. 60). 
For free space we shall have g=0, and the following numerical values of the 
fundamental constants (p. 82): 
impedance of free space mo ~ 1207 ohms, (9-4) 
characteristic velocity vo ~ 3.108 meters/second. 
a * * ao °K *” x co 

Surprising as it may appear to transmission engineers and sound engineers, 
who daily handle their respective characteristic impedances Zo or pc, there 
still are very competent physicists who balk at the idea of free space having a 
characteristic impedance of about 377 ohms. Yet, in the words of Professor 
Ronold W. P. King:! “The existence of such a characteristic resistance for 
electromagnetic effects is just as mysterious, but not more so, than the exist- 
ence of the finite velocity vo.” Dr. Schelkunoff explains very well how this 
constant could have been overlooked by the builders of the classical theory: 
“The physicist concentrates his attention on one particular wave: a wave of 
force or a wave of velocity or a wave of displacement. His original differential 
equations may be of the first order and may involve both force and velocity; 
but by tradition he eliminates one of these variables, obtains a second order 
differential equation in the other and calls it the “wave equation.” Thus he 
loses sight of the interdependence of force and velocity waves... ” (p. vii). 
Still, it is surprising to see that one has started with two constants €9 and yo, 
recognize the fundamental importance of their product, and not enquire about 
their ratio. 

Then, the reader will ask, how can the Theory of Relativity give a leading 
role to the velocity of light and not mention the impedance of free space. 
Has Einstein no use for ny? Well, he has, and he has not. First, an essential 
point in Special Relativity is the merging of the magnetic and the electric 





1 Mimeographed “Notes on Antennas” for the course of Electronics and Cathode Ray 
Tubes (Eng. 270), Harvard University. 
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fields into one skew-symmetrical tensor. When doing this in the MKS sys- 
tem, homogeneity requires the use of the components of E and of nof7; but 
the factor mo is not apparent, for instance, in the equations on p. 44 of “The 
Meaning of Relativity” (by A. Einstein, Princeton Univ. Press, 1923) which 
uses a system of units in which no=1. Secondly, if we try to connect the uni- 
versal constant 7» with other members of this interesting family, we find that 
no times a (charge)? has the dimensions of “action,” and more precisely that 


2h 
137 


(e=charge of the electron, k= Planck’s constant). We see from this that there 
is more to no than appears in Special Relativity, the first step in the successive 
Einsteinian extensions of Maxwell's theory. 


noe? = 


* * * *~ * * * * ; 

We have dealt at length with this question of the “impedance of free space” 
because it exemplifies the spirit of the whole work. It occurs in the course 
of a short but apt presentation of the “Fundamental Electromagnetic Equa- 
tions” (Chapter IV), immediately applied to harmonic oscillations. The book 
as a whole is devoted not to Electromagnetism in general but, as specified in 
the title, to Electromagnetic Waves. 

Three preliminary chapters introduce the more advanced mathematical 
tools which will be used, but sparingly, in what follows: such topics as contour 
integration, Bessel and Legendre functions. Chapter V is a short and original 
presentation of Network Theory. 

The central part of the book begins with Chapter VI, “About Waves in 
General,” a sort of preview of the questions which will be treated in detail 
later, during which we are introduced to radiation from given currents, 
propagation along wave guides, and to such general tools as electric and mag- 
netic current sheets, the method of images and conformal representation. 

In the following four chapters, we meet the most thorough treatment 
available of the propagation of waves, guided or bounded, in one, two and 
three dimensions. It is impossible to do justice here to the richness of the 
material, which must have cost tremendous labor and which is in great part 
taken from the author’s own publications. We find in Chapter IX, however, 
classical problems of Fresnel optics, adroitly adapted to contemporary radio 
needs. Chapter XI is a relatively short treatment of antenna theory, princi- 
pally of conical antennas, and in the last chapter we return to wave guides 
and solve various problems involving discontinuities, even to an iris or a 
transversal wire. This subject is still under development by the author, and 
the readers of the Quarterly have had the benefit of one of its recent exten- 


sions.? 


2 See Quarterly of Applied Mathematics, 1, 78 (1943). 
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The specialist in wave propagation has no need to be told of the value of 
this book; but the reviewer would like to explain to his fellow non-specialists 
why it is particularly important that they should not miss it. When the re- 
sults of much present-day research will suddenly be made available, it will 
be a hard task to catch up, not only with the new knowledge, but still more 
with the new modes of attack. The borderland between radio and optics is 
one of the fields from which great things can confidently be expected. Dr. 
Schelkunoff’s book is a great opportunity for those not at present engaged in 
research to get familiar with methods which they will want to use tomorrow. 

P. LE CORBEILLER 


























